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ABSTRACT 


To reduce the number of phase sampling comparators in a direction finding (DF) 
interferometer antenna, a new array based on a robust symmetrical number system 
(RSNS) is described. The RSNS is used to decompose the spatial filtering operation into 
a number of parallel sub-operations that are of smaller computational complexity. Each 
sub-operation (interferometer) symmetrically folds the phase with folding period equal to 
2Nm\ where N is the number of channels that are used and m\ is the channel modulus. A 
small comparator ladder mid-level quantizes each folded phase response. Each sub¬ 
operation only requires a precision in accordance with that modulus. A much higher DF 
resolution is achieved after the N different RSNS moduli are used and the results of these 
low precision sub-operations are recombined. The parallel use of phase waveforms 
increases the antenna resolution without increasing the folding rate of the system. The 
new antenna is constructed and tested in an anechoic chamber, and the results are 
compared with the experimental results of a previously tested optimum symmetrical 
number system (OSNS) array. Although the dynamic range of the RSNS is somewhat 
less than the OSNS, the inherent Gray code properties make it attractive for error control 
in phase sampling interferometry. 
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I. INTRODUCTION 


A. PHASE SAMPLED INTERFEROMETRY 

Direction finding (DF) systems find applications in many disciplines including 
geo-location, navigation, and targeting. Law enforcement officials, military personnel 
and wildlife managers use DF systems extensively for tracking. Radio DF systems 
commonly use phase sampled linear interferometer arrays to measure the direction of 

i 

arrival of a single incident planar waveform in order to estimate the geographical bearing 
of the emitter. Each interferometer consists of two receive elements separated by a 
distance d (baseline). A waveform impinging on an ideal interferometer has the same 
amplitude at both element locations. Thus all of the information about the direction of 
arrival is contained in the phase of the signal measured by the different array elements. 
Consequently, most existing radio DF systems are based on phase-only measurements. 
As the planar waveform impinges upon the receive elements, the difference in phase is 
measured, and the direction of arrival is estimated. Ambiguities in the phase difference 
occur if the separation between interferometer elements is greater than half the electrical 
wavelength. These ambiguities cause ambiguous estimations of the direction of arrival. 

To resolve the ambiguities, the addition of a second interferometer (third element) 
can be used [1,2]. Multiple baseline systems can become quite long when high resolution 
is required. This long baseline may be unacceptable in volume confined applications 
such as spacecraft, unmanned aerial vehicles, and hand held DF systems. Many 
alternative approaches to resolving the ambiguities within a phase sampling system have 
been recently investigated. 

The extended phase interferometry technique [3,4] incorporates both calibrated 
phase and amplitude response data from the antenna arrays. This technique appropriately 
weights the square of the baseline phase differences with the antenna gains. The 
incorporation of amplitude data provides significant performance improvement over 
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phase-only interferometry; however, it requires a modest increase in computational 
complexity. 

Pseudo Doppler techniques [5] have also been investigated in the 2-2000 MHz 
range to prevent ambiguities. This system employs a circular array of four omni¬ 
directional antennas. A commutation switch samples the antennas such that the resulting 
phase modulation on the received signal encodes the direction of arrival. As the elements 
are commutated a FM discriminator responds with impulses with size approximately 
proportional to the phase step between each element. The maximum spacing between 
adjacent elements is less than half a wavelength. The phase of the output is used to 
estimate the angle of arrival. The single channel interferometer minimizes cost and 
power consumption by avoiding a multitude of expensive receivers and by employing the 
cost-effective pseudo Doppler antenna. 

The symmetric pair antenna array [6] provides accurate DF by means of low- 
resolution phase measurements, achieved by virtue of the symmetry of the array 
geometry. This type of array consists of pairs of antenna elements that are arranged 
around a common phase center. Due to symmetry, any arriving wavefront generates a 
phase lead at one member and a phase lag of equal magnitude at the other member of the 
symmetric pair. Since the polarity of the output signal is dependent on the angle of 
incidence, it can be used to determine the signal direction of arrival. The sign bit pattern 
of the pairs of a given symmetric pair array will change for different angles of incidence 
and, assuming an unambiguous arrangement, (annular zones on the surface of a sphere) 
will be distinguishable. 

This thesis presents a phase sampled linear interferometer array that incorporates 
a robust symmetrical number system encoding. The system has no moving parts or 
switches. Due to the inherent Gray code (adjacent codes differ by one bit) properties, the 
robust symmetrical number system (RSNS) efficiently resolves the ambiguities and 
resists mutual coupling induced errors. RSNS is used to decompose the spatial filtering 
operation into a number of parallel sub-operations that are of smaller computational 
complexity. Each sub-operation (interferometer) symmetrically folds the phase with 
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folding period equal to 2Nm { for the Robust Symmetrical Number System (RSNS) where 
N is the number of channels that are used and m t is the channel modulus. A small 
comparator ladder mid-level quantizes each folded phase response. Each sub-operation 
only requires a precision in accordance with that modulus. The low precision demanded 
in each sub-operation permits the system to overcome phase errors induced by mutual 
coupling. A much higher DF resolution is achieved after the N different RSNS moduli 
are used and the results of these low precision sub-operations are recombined. The 
parallel use of phase waveforms increases the antenna resolution without increasing the 
folding rate of the system. Previous experimental results from an Optimum Symmetrical 
Number System (OSNS) prototype antenna revealed encoding errors at the code 
transition points. Although the dynamic range of the RSNS is somewhat smaller than 
OSNS, the inherent Gray code properties of the RSNS make it attractive for error control. 
A prototype RSNS system has been designed and built to investigate the effectiveness of 
this signal processing technique. 

B. PRINCIPAL CONTRIBUTIONS 

The focus of this thesis research has been to design, fabricate, and test prototype 
antenna arrays that investigate the performance of the RSNS architecture. This antenna is 
the first of its kind incorporating innovative microwave technology (stripline antenna 
elements), developing new digital hardware designs, and deriving unique performance 
equations. The antenna design results in a system that is resistant to phase errors while 
simultaneously providing high resolution with a minimal amount of hardware. 

The design efforts were initially focused upon the phase detection process. 
MATLAB programs were written to simulate the design. Several array configurations 
were built and evaluated. MATLAB programs were also written to evaluate the 
experimental results. The decision to build an array with a 6-bit resolution restricted the 
element spacings and demanded the design and fabrication of innovative stripline dipole 
antenna elements and associated microwave circuits. 

The sub-operations consist of mid-level quantization accomplished by parallel 
comparator networks. The derivation of the comparator thresholds enabled the proper 
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quantization to occur. The comparator networks process the folding waveform output of 
the phase detection circuit. The proper phase relationship between the antenna elements 
was derived and verified through MATLAB simulation. 

The initial results from the array possessed phase errors attributed to mutual 
coupling between the antenna elements. A scale factor was derived and performance 
simulated in order to increase the element spacing and reduce the mutual coupling 
effects. 

A hybrid of Gray code and binary EEPROM mapping enabled optimization of the 
digital circuit design and reduced the required number of components. A MATLAB 
program simulated the performance of the digital circuit. 

Data processing was conducted by collecting the EEPROM output and phase 
detector outputs in parallel with a LAB VIEW program. Over 100 patterns studied the 
impact of mutual coupling, antenna scaling, amplifier harmonics, and signal multi-paths. 

C. THESIS OUTLINE 

In Chapter II, two-element interferometry theory is reviewed. The equations 
governing the sampling of an incident wave front are presented and the variables that 
effect the output are discussed. 

In Chapter III, the Optimum Symmetrical Number Systems is presented. The OSNS 
prototype antenna is described, and the experimental results are presented. 

In Chapter IV the application of the RSNS theory to phase sampled interferometry is 
described. 

Chapter V derives the design equations and application to the specific prototype 
antenna. The geometric derivations of the required phase and comparator settings are 
shown. The MATLAB simulations of the prototype antenna arrays are provided. The 
folding waveforms and transfer functions used to design the prototype are presented. 

Chapter VI describes the hardware design and anechoic chamber testing procedures. 

Chapter VII presents the experimental results from the prototype antenna arrays. 
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Chapter VIII summarizes the results and conclusions derived from the RSNS 
prototype antenna array. It also presents areas for further research and insights into the 
possible results of that research. 

Appendix A presents the MATLAB codes used to design, simulate, and evaluate the 
RSNS antenna. 

Appendix B shows the LABVIEW codes used for data collection. 

Appendix C describes the array design, predicted folding waveforms, and predicted 
transfer function of a 3-channel RSNS antenna used to investigate the folding waveform 
phase errors caused by mutual coupling. 
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II. PHASE SAMPLED INTERFEROMETRY 


A. DIRECTION FINDING SYSTEMS 

Radio direction finding is a class of direction finding by which the direction to a 
radio source is determined by means of a radio direction finder (RDF). A RDF is a 
passive receiving system that measures the energy extracted from the passing 
electromagnetic radio wave to obtain direction-of-arrival (DOA) information. 

A radio direction finding system consists of four essential functional elements [7]. 
The antenna extracts the electromagnetic energy and converts it to a signal containing 
DOA information. The receiver down-converts, amplifies, and processes the antenna 
signal to intermediate frequencies (IF) or baseband. The post-receiver processor further 
processes the signal to obtain the DOA information. The DF information processing- 
read-out-display unit prepares the basic DOA data for transmission to users of the DF 
information. The antenna is the key subsystem in a RDF. It extracts energy from the 
incident electromagnetic field and provides output signals containing the incident energy 
DOA information [8], 

At microwave frequencies, antennas typically have aperture sizes of several 
wavelengths. Small aperture RDFs determine direction-of-arrival information by one of 
the following three measurement methods [7]: 

1. Amplitude response, 

2. Phase difference, 

3. Time delay. 

More specifically, small aperture DF systems may be separated into four functional 
categories based on the fundamental measurement methods presented above. The four 
basic categories are as follows [7]: 

Category I: Systems using either direct or comparative amplitude response of the 
antenna subsystem for DOA information. 
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Category II: Systems using the phase difference between disposed antenna 
elements with the phase differential converted to amplitude DOA information. 

Category III: Systems using the phase difference between disposed antenna 
elements for DOA information. 

Category IV: Systems using the time-of-arrival difference between disposed 
antenna elements for DOA information. 

The RSNS antenna investigated in this thesis falls into Category III. 


B. PHASE COMPARISON SYSTEM 

A two-element linear interferometer is shown in Figure 2-1. The two antenna 
elements are spaced a distance d apart, and the incident plane wave arrives with bearing 
angle 6. In this phase monopulse configuration the angle 6 is measured from the 
perpendicular to the baseline axis and can take on values n H> 9 > -n 12. The 
electromagnetic wave arrives first at antenna element 1 located at x = d/2. Without loss 
of generality, the phase of the wave arriving at the first element at time t is , = Inft . 
The plane wave arrives at antenna element 2 (located at x=-d/2 ) after traveling an 
additional distance of d sin(#). The phase of the wave arriving at the second element is 
then 


2 = 2 xft-—dsm{9). (2.1) 

The phase difference between the two elements is 

2n 

= —dsm{9 B ) = kd sin(0), (2.2) 

A 


where k is the wave number, is a function of the incidence angle of the wave. 
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The signals received by the antenna elements are phase detected (multiplied 
together and processed by a lowpass filter) resulting in an output signal whose frequency 
is the difference between the two input signal frequencies. Because the signals have the 
same frequency, the mixer output is a DC voltage whose value is proportional to the 
difference between signal phases y/ x and y / 2 . This difference is not purely Ay/ because 
of the time delays due to the different transmission line lengths from each antenna 
element. However, these line lengths are known and defined as <f> x and <p 2 . The phase 
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error due to the line length differences is 0 = 0[ —(f> 2 and can be compensated for in the 
angle estimate. Let the signals from the two antenna elements be 

v i (0 = V cos[^ + 0 t (/)] (2.3) 

and 

v 2 (0 = V cos[^ 2 + 0 2 (03, (2.4) 

where V is the maximum value of the voltage at the antenna elements. The lowpass 
mixer output voltage is 

V 2 y 2 

v ou< = —cos(A If/ + 0) = —cos (kd sin(0) + 0) (2.5) 

which contains the plane wave angle of arrival (AOA) information. For values of 
d = X!2, 0 = 0, and A^ = ^sin(<9), as the AOA 6 varies from -n/2 to nil , the 
phase difference A0 varies from -n to n as shown in Figure 2 - 2. The output voltage 
from the phase detector is also a function of the phase difference and is a symmetrical 
folding waveform. Together these relationships give the phase detector output voltage as 
a symmetrical function of the AOA as shown in Figure 2 - 3 for ^ = 2/2. 

Ambiguities are generated for baselines where d>XI 2. That is, the phase 
detector output voltage is highly ambiguous with a single phase corresponding to many 
angles of arrival. The number of folding periods n that occur within an AOA of n 
radians is 


2d 

T' (2-6) 
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Angle of arrival (degrees) 



Figure 2-3: Mixer Output Voltage Versus Angle of Incidence, d= X /2. 


11 










For example, with d = 1.5A>, n =15 folds are available as shown in Figure 2-4. 
The folding period is not constant but grows larger in proportion to the off-perpendicular 
angle 0 because of the sin(<9) dependence in (2.5). 



Figure 2-4: Mixer Output Voltage Versus Angle of Arrival for d =7.52. 


The ambiguities within the symmetrical folding waveforms represent the phase 
difference between the elements and can be resolved by using additional interferometers 
in the linear array. Typically, each interferometer in the linear array symmetrically folds 
the phase response with the folding period between interferometers being a successive 
factor of two, or d 4 = 2</ 3 = 4 d 2 = 8 d x . High-speed binary comparators are used to 

produce a digital output. The folding waveforms are shifted appropriately using a phase 
shifter in each channel to achieve a Gray code result. The folded output from each phase 
detector is then quantized with a single comparator with a normalized threshold level 
T=0.5. Together, the comparator outputs directly encode the signal's AOA in the Gray 
code format. This approach makes use of the periodic dependence of the interferometer’s 
phase response on the applied plane wave's AOA and the distance between the elements 
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of each interferometer. One of the major limitations associated with this approach is the 
achievable resolution. For each folding period to be a factor of two longer than its 
predecessor, the distance between the elements must also be doubled. That is, an 8-bit 
DF antenna using the previous scheme would require element spacings 
A/4,A/2,A,2A,4A...,32A with a total baseline length of 32A . This distance doubling of 

the element spacings requires complex analog hardware, adversely affects the physical 
implementation of the DF architecture, and ultimately constrains the achievable 
resolution. 
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in. OPTIMUM SYMMETRICAL NUMBER SYSTEM 


DIRECTION FINDING ANTENNAS 

A. OPTIMUM SYMMETRICAL NUMBER SYSTEM 

Incorporation of the optimum symmetrical number system (OSNS) encoding 
overcomes the major limitations of the binary encoding approach and provides an 
efficient method to enhance the resolution of the array while minimizing the total 
hardware baseline. 

The OSNS is composed of a number of pairwise relatively prime (PRP) moduli 
m r The integers within each OSNS modulus are representative of a symmetrically 
folded waveform with the period of the waveform equal to twice the modulus, i.e., 2 m i 
[9]. For m given, the integer values within twice the individual modulus are defined by 
the row vector 


x m = [o,l,---,m-l,m-l,---,l,o]. (3.1) 

Figure 3-1 shows part of the OSNS folding waveforms and the integer values within the 
modulus for m\=5 and rri 2 = 6 . The horizontal axis represents the normalized input. The 
vertical lines represent folding waveform reference levels. The numbers in squares at the 
top of the figure represent the number of reference levels that are crossed by the folding 
waveform for a given input value and correspond to the row vector x m . The period of 

one complete fold is equal to 2m, =10 and 2m 2 = 12. 

Due to the presence of ambiguities, the integers within the vector (3.1) do not 
form a complete system of length 2m by themselves. The ambiguities that arise within 
the modulus are resolved by considering the paired values from all channels together. If 
N pairwise relatively prime moduli are used, the dynamic range (the number of 
unambiguous vectors) of this scheme is [10]: 
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(3.2) 


N 


M = Um i . 

i=i 


T = I 
T = 2 
T = 3 


T = 4 


T ■ I 
T = 2 
T = 3 
T = 4 
T = 5 



Figure 3 - 1: Symmetrical Folding Waveform for OSNS m t =5 and m 2 =6. 


m 2 = 6 


This dynamic range is also the position of the first repetitive moduli vector. As shown in 
Figure 3 - 1, the dynamic range for them, =5 and m 2 = 6case is M= 30. That is, no set 
appears twice for a normalized input in the range [0,29]. 

B. OSNS ANTENNA ARCHITECTURE 

Figure 3-2 shows the schematic diagram of a N=2 channel OSNS DF antenna. 
The array consists of three elements, with both channels sharing a common element. 
With the correct distances between elements, each interferometer folds the phase 
response at 2m,. The phase response from each channel is quantized using m,-l 
comparators. The integers shown at each step on Figure 3 - 1 represent the number of 
comparators ON due to the phase voltage exceeding the comparator matching threshold 
voltage. 
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Figure 3-2: Optimum Symmetrical Number System Antenna Architecture. 


The distance between each pair of elements is a function of the modulus, m,. The 
phase response must have the correct number of folds to cover the dynamic range. For a 
given modulus the number of folds within the dynamic range is 


n- t = M / 2 m j . 


(3.3) 


Since each fold of the array corresponds to a change in phase of 7t for d = /i/2, 
the required spacing between the reference element and the respective element of channel 
i is given by 









The bin widths are uniform in size in sin (0) space with a width 2/M. The spatial 
resolution will be derived for Equation (4.18) and is given by: 


= arcsin 

v 


2k-M+ 2 
M 


\ 

J 


-arcsin 


2 k-M 
M 


\ 

J 


(3.5) 


where k e {0,1,--• A/ — l} and k = Ois the bin that starts at 0 = -n/2. Note that the spatial 
resolution is not uniform over the range of -n/2 to n/2. The resolution is finest at 
broadside (0 = 0 °) and increases with angle off broadside, because of the 
sin(0) dependence in (2.5). 

In some applications element placement at certain locations on a platform may be 
prohibited, in which case an alternate OSNS design could be used to provide the same 
unambiguous AOA resolution, using e.g., a different number of comparators, as shown in 
Table 3-1. When the channels are recombined the result is a high accuracy 
unambiguous estimate of the AOA. 


OSNS Moduli 

Element Spacings 

Comparators Required 

mi = 3, m 2 = 20 

dj = 5X, d 2 = 0.75X 

21 

mi = 4, m 2 = 15 

di = 3.75X,d 2 = IX 

17 

mi = 5, m 2 = 12 


15 

mi = 3, m 2 = 4 

di = 5X, d 2 = 3.75X 

9 

m 3 = 5 

d 3 = 3X 



Table 3 - 1: Several OSNS Arrays Possessing the Same Spatial Resolution (M = 60). 


C. EXPERIMENTAL RESULTS 

An OSNS array based on the moduli m\ = 6 and m 2 = 11 was designed, fabricated 
and tested at a frequency 8.5 GHz [11]. The radiating elements are open-ended 
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waveguides as shown in Figure 3-2. For moduli m\ = 6 and m 2 =11, d\=2>.Z2 inches 
and ^2=2.08 inches. With M=66, the number of folding periods n\=5.5 and «2 = 3. To 
provide an adequate signal-to-noise ratio, a low-noise amplifier is included at the output 
of each interferometer element. Since the common element splits the signal into N paths, 
an attenuator is placed in the other branches to balance the amplitudes. The amplifiers 
operate in saturation so that the mixer input signal levels are independent of the angle of 
incidence. A fixed phase shifter is also included in one branch of each interferometer so 
that the symmetrically folded phase response waveforms from each mixer may be 
aligned. This alignment insures that the comparators in the digital processor properly 
sample the phase waveform and encode it in the OSNS. 

The anechoic chamber facility at the Naval Postgraduate School was used for the 
antenna pattern measurements for each pair of elements. The measured mixer outputs are 
shown with the ideal folding waveforms in Figure 3-3 and Figure 3-4. Since the array 
is not useful in the endfire regions, the phase resppnse from each channel was minimized 
at input AOA of -50°. This is the point where the digital code starts. Introducing the 
appropriate phase shifts to the signals at a point between the elements and the mixer 
controls the location of the voltage minima. A simulation of the transfer function for the 
OSNS array is shown in Figure 3-5. The measured transfer function (output of digital 
hardware) is shown in Figure 3-6. 

The spikes in the transfer function are AOA encoding errors. The phase 
waveforms must cross the comparator thresholds (code transitions) at the same time 
across all channels. The OSNS can present a problem at each of these specific AOAs. 
When the set of comparators does not change state synchronously, the recovered 
amplitude will have a large error. Encoding errors can also occur due to imperfections in 
the transmission line phases. These errors however, can be isolated with additional 
comparators and corrected using additional circuitry as described in [9]. An alternate 
approach however is to develop a symmetrical number system that eliminates the 
possibility of this type of encoding error. That is, a symmetrical number system where 
only one comparator threshold is crossed at any one particular code transition. 
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Mixer Output Voltage (Volts) 



Figure 3-3: Comparison of Measured and Predicted Modulus 6 Mixer Output [from 10] 



Figure 3 - 4: Comparison of Measured and Predicted Modulus 11 Mixer Output [from 10] 
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Figure 3-5: Simulated Transfer Function of OSNS Measured Mixer Output [from 10,11], 


Experimental Results 



Figure 3-6: Measured Transfer Function of OSNS Prototype Antenna [from 10,11] 
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IV. ROBUST SYMMETRICAL NUMBER SYSTEM 
DIRECTION FINDING ANTENNAS 


A. ROBUST SYMMETRICAL NUMBER SYSTEM 

In the robust symmetrical number system (RSNS), N different periodic 
symmetrical waveforms are used with pairwise relatively prime integers mi, m 2 , ..., m^. 
The RSNS is based on the following sequence 

Xmi “[0,1,2, .... m- 1 , m, m- 1 ,..., 2,1] , (4.1) 

where x m i is a row vector and m is a positive integer (m > 0). In an iV-channel RSNS, 
where N >2, the basic sequence for the i th channel (modulus m) is 

x m 2=[0, 0,..., 0, 0, 1, 1,..., 1, 1,..., m, m,..., m, m, ..., 1, 1,..., 1, 1J. (4.2) 

In this sequence each value in the x m i row vector is entered N times in succession. This 
sequence is repeated in both directions, forming a periodic sequence with the period 

Prsns = 2m iN. (4.3) 

Considering a single channel, the discrete states of the robust symmetrical number system 
can be expressed as [12] 




2Nm t +N-n + Sj—l 
- ^ 


Sj <n< Nm , + s, +1 
Nm i + Sj + 2 < n < 2 Nm i + 5 , -1 


(4.4) 
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where g is the n th term of channel i, m l is the channel modulus, s t is a corresponding 
sequence shift s,- = 0, 1, 2,..., N-l (mod N) and N > 2 is the number of channels in the 
system. The values (sj, S 2 ,..., sn} must form a complete residue system modulo N. 
Because of the relative property of the shifts, one of the shift values will be set equal to 
zero. The index n corresponds to the input value. The discrete states of the RSNS are 
indexed using the s i = 0 row vector with the index starting from the first zero. Note that 
the largest integer within each periodic sequence is the modulus m. 

An A-channel RSNS is formed of vectors by picking N moduli m h and N shift 
values Si, 1 <i <N. Since the fundamental period for channel i is 2Nm u the period for the 
RSNS vectors must be a multiple of 2Nmj. Therefore the fundamental period for the 
RSNS is 


PFrsns ~ [2miN, 2 m 2 N,2m^N], 


(4.5) 


where [a } , a 2 , .... a N ] is the least common multiple of a h a 2 , .... a N . From number theory 
[13], 


PFrsns = 2N[mi, m 2 ,..., m N J. (4.6) 

The system dynamic range (m rsns ) of the RSNS is the maximum number of 

distinct vectors without a redundancy. For two channels, with the moduli separated by 1 
or 2 [14], 

M RSNS = 4m i +2m 2 -5. (4.7) 

For two channels, with the moduli separated by 3 or more, 

Mrsns = 4m i +2 m 2 ~2. (4.8) 
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Checking against computer results when mi=8 and W2=17, M = 64 which fits (4.8). 

B. RSNS ANTENNA DESIGN 

Element spacing is a fundamental characteristic of any antenna. The RSNS 
architecture is capable of high resolution with very small baselines. The choice of 
spacing is not arbitrary; it defines the folding waveforms produced by the analog circuit. 
The element spacing is a function of the modulus, number of channels, dynamic range, 
and signal frequency. The spacing is found in three steps: determine the modulus folding 
period, determine the number of folds contained in the dynamic range, and factor in the 
frequency. 
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Table 4 -1: Vector Series for Two Channel RSNS, Modulus 3 and 4. 


Table 4-1 shows the integer values within the RSNS for mi=3 and m2=4. Recall 
from (4.3) that the folding period depends on the modulus m, and the number of channels 
N and given by 

Prsns = 2m,iN. (4.9) 


The number of folds is: 

A A 

M _ M 
Prsns 2m t N 


(4.10) 


The folds must fill the antenna’s field of view (FOV). For the RSNS antenna, the FOV is 
180°. In phase space, 180° corresponds to XU. Therefore, the required spacing, d, is 
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(4.11) 


, X MX 

a, =n, — = -. 

2 4m,N 

Comparing (4.11) to (3.4), for systems with the same dynamic range, the RSNS element 
spacing is smaller than the OSNS element spacing by a factor of N. 

C. COMPARATOR THRESHOLDS 

In the RSNS, m { comparators are required for each channel. Consider the cosine 
of 6 for -n < 9 < 0, shown Figure 4-1. This represents a half-fold of a folding 
waveform. If the folding waveform represents the modulus 3 channel of a two-channel 
system, the modulus 3 component of the RSNS vector string is {... 1,0,0,1,1,2,2,3,3,2...}. 
The first zero and second three belong to the adjacent half-folds. This leaves six bins for 
this half-fold. To space the bins equally, each bin must span tc/6 radians. The modulus 
3 component changes when a comparator changes state. The three comparator thresholds 
for this case are cos(-tt + f), cos(-?r + ^), and cos(-tt + f). Figure 4 - 2 shows these 
comparator states. 



Figure 4 -1: Half Fold of Cosine Theta. 
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Figure 4-2: Modulus 3 Comparator Thresholds. 

The expression for the y'-th normalized comparator threshold in a modulus tm, 
channel is observed graphically to be: 


V = cos 

J 



(4.12) 


The comparator thresholds depend only upon the channel modulus. This powerful 
relationship is required prior to determining the proper phase relationships of the folding 
waveforms. 


D. CHANNEL PHASING 

A computer search algorithm determines the dynamic range of a set of moduli. 
Equations exist for special cases such as Equations 4.7 and 4.8. Even when the value of 
the dynamic range is known, the computer search algorithm is used to determine the 
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location of the dynamic range within the vector string [14]. For the case of mi =3 and 
w 2=4, the dynamic range is 15 and falls between vectors 6 and 20 as shown in Table 4 - 
2. The bin number is a re-indexing of the vector, counting the first unambiguous vector 
as zero. 
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Table 4-2: Dynamic Range of RSNS System, Modulus 3 and 4. 


The OSNS has folding waveforms that are either even or odd functions about 
broadside. For RSNS, the first vector is not the zero vector, and the dynamic range is not 
a multiple of the fundamental period. Coincidental cases of symmetry may occur in the 
RSNS dynamic range. This is clearly seen in Table 4-1, which has the dynamic range 
shaded. The ambiguous vectors in vectors 5, 9 and 21 delimit the dynamic range. The 
first vector of the dynamic range is vector 6, is (3,3). Neither modulus has symmetry 
about the midpoint, vector 13. The lack of symmetry in the RSNS dynamic range makes 
it necessary to determine this unique phase value at broadside for each channel. From 
(2.5) the normalized phase detector output is 

mx(6) = cos(kd s,mO + <j>). (4.13) 

When 0=0, sin (9=0, so mx(#=0)=cos(^). Theta equal zero corresponds to the center of 
Bin 7. This is broadside incidence. At the left edge of Bin 7, the folding waveform is at 
its minimum, and at the right edge of Bin 10, the folding waveform is at maxim um This 
is exactly the waveform of Figure 4-2. 

Expanding the lower left comer of this graph, Figure 4-3 shows bins 7 and 8. 
The left edge of Bin 7 corresponds to <jf =-n and the right edge of Bin 7 is the comparator 
threshold, <j) = =^-. The center of Bin 7 is midway between these angles: <j> % - . The 
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sign is negative because the phase is on the negative half cycle. A similar phase 
determination for Modulus 4 shows the left edge of Bin 7 is at the comparator #2 
threshold, which is (from Equation 4.11) (j> = ^-. The right edge of Bin 8 is at the 

comparator #1 threshold, = The center of Bin 7 is one quarter of the angular 
distance between the thresholds, or <f> 4 = ^. The signs of these angles is positive 

because the folding waveform is in the decreasing half cycle. This phase difference 
between the elements can be inserted with a phase adjuster in the analog microwave 
circuit. The resulting folding waveforms with the correctly determined phases is shown 
Figure 4-4. The equations that describe the folding waveforms are: 


f 

mx 3 (0) = cos(kd 3 sin0 + <f> 3 ) - cos ■ 

V 

f 

mx 4 (0) = cos(M 4 sin# + ^ 4 ) = cos 
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Note that (4.14) and (4.15) are independent of frequency. 


(4.14) 

(4.15) 
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E. MAPPING RSNS DIGITAL OUTPUT TO ANGLE OF ARRIVAL 


Whenever a folding waveform exceeds a comparator threshold, the comparator 
turns ON. For the folding waveforms shown in Figure 4 - 4, the numbers of comparators 
in the ON state are shown in 

Figure 4-5 (integer values within RSNS code). The active comparator map 
exactly matches the RSNS unambiguous vector string shown in Table 4-2. The 
dynamic range of this string falls in vectors 6-20. This is the string of vectors that are 
useful for DF, so those bins are indexed as shown in the bottom row of Table 4-2. 

Re-indexing the vectors in the dynamic range with a bin number enables the 
RSNS code to map the AO A. The comparator states of 

Figure 4 - 5 are unambiguously mapped into a bin as shown in Figure 4-6. 



1 -80 -60 -40 -20 6 .20 40 60 80 

Angle of Arrival (degrees) 

Figure 4-5: Active Comparators. 
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Figure 4-6: Bin Number to AOA Mapping. 

The folding waveforms of Figure 4 - 4 expand as a function of 6 . This causes the 
bin widths (the angle encompassed within a bin) in Figure 4 - 6 to expand as |<9| increases 

(the incident wave diverges from broadside). This effect is caused by the embedded 
trigonometric functions of Equations 4.14 and 4.15. This effect is obvious when 
comparing the angle included in Bins 7 and 14. The bin widths are expanding as 1/sin 9 , 
as described in Section F. In Figure 4-7, the bin mapping is shown as a function of the 
sine of the angle of arrival. In this graph, the bin widths are uniform across the sin# 
domain. 


32 

















Bin Number 






For the first bin (Bin 0), the center is located as sin<9 0 =-l + BW/2, and the center of the 
last bin (Bin 14) is located at sin# 14 = 1-BW/2 . A consequence of quantizing the angle 
of arrival into a bin is that the RSNS system reports any signal that falls within a bin as if 
it arrived at the bin center, 9 k . That is to say, for one angle within each bin, the 

estimation is exact, but for the remaining angles a reporting error exists. Figure 4 - 9 
illustrates the cause of this sampling quantization error. The quantization error will be 
analyzed in the next chapter. The bin width in 9-space is 


r. = arcsm 

k 


/ A A \ 

zfc-M + 2' 

A 

M 


- arcsm 


/a A \ 

r 2k-M 




M 


(4.18) 


The best resolution is at broadside (6* = 0) and becomes worse in the end-fire 
(j#| —> 7t 12) regions. Taking the arcsine of each side of Equation 4.17 gives the 
quantized angle of arrival 9 as: 


9 ; = sin 

k 


-1 


/ * 

Z& + 1 

A 

M 


-1 


(4.19) 


This is the RSNS transfer function. 
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Figure 4-8: Reported Angle of Arrival. 



Figure 4-9: Quantization Error for Three Bit Resolution. 
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F. 


SCALE FACTOR AND MAPPABLE FIELD OF VIEW 


The number of folds «, is the number of folds the a given modulus with a 
specified dynamic range is capable of mapping. As the element spacing d t is increased, 
the number of folding waveforms increases. Increasing d compresses the existing 
waveform and symmetrically adds additional folds at the extreme angles. This 
compresses the aperture that the RSNS maps. The maximum mappable aperture, 0 mm , 
that the RSNS dynamic range is capable of mapping is determined by a scale factor, £ 

The number of folds is a linear function of d (Equation 4.11). For d=2/2, the 
mixer output has one fold. The ratio of the element spacings where d is the unsealed 
separation and d' is the scaled element separation is the scale factor. 


d_ _ sin# 
~d’ ~ sin#' 


(4.20) 


The mixer output function of 0 is shown in Figure 4-10 for three scale factors. The mixer 
output as a function of sin# is shown in Figure 4-11 for the same three scale factors. 
Including the scale factor in the element spacing, it follows that 


m 

ANmg ’ 


(4.21) 


The new angle resulting from the transformation is: 


O' = arcsin 



(4.22) 
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Normalized Mixer Output Normalized Phase Output 



Figure 4 -10: Effect of Three Scaling Factors as a Function of Theta. 



Figure 4-11: Effect of Three Scaling Factors as a Function of Sine Theta. 
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For example, for 6 = nil and -nil, the new angles O' are shown in Figure 4-12. 
These particular angles represent 0 mm the maximum mappable aperture of the scaled 
waveform. If the separation of all elements is scaled by the same scale factor, all of the 
resulting folding waveforms are compressed similarly. The RSNS mapping is 
maintained; however, the mapping is contained within 0 mm . The scaled antenna 
determines the angle of arrival according to: 


f 


6 - arcsin 


V 


lk + l 
M4 



(4.23) 


Beyond 0 mm , the folding waveforms continue to be processed by the RSNS circuitry, but 
the resulting vectors either do not map into the dynamic range or are ambiguous as shown 
in Figure 4-13. For comparison, the same transfer function operating on an unsealed 
antenna produces the transfer function shown in Figure 4-14. 



Figure 4 - 12: Determination of Maximum Mappable Aperture. 
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There are several factors affected by changing the spacing of the elements. The 
increased spacing reduces the antenna’s useful field of view, but the larger spaces can 
decrease mutual coupling. The new field of view reduces the mappable aperture, but it 
also improves the antenna resolution. The folding waveforms at the endfire do not 
behave properly and are discarded in the scaled antenna. 

The process described in this chapter demonstrates the method of determining the 
transfer function for a given RSNS antenna configuration. In the next chapter these 
methods are used to design a high resolution antenna. 
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V. PROTOTYPE DESIGN AND PREDICTED RESULTS 

A. DETERMINATION OF MODULUS AND DYNAMIC RANGE 

The previous chapter presented the design equations for an RSNS architecture. This 
chapter will apply those equations to the design of a prototype antenna array. Two arrays 
are presented. The arrays are identical except that one is an unsealed array (£ = 1) and 

the other is scaled with £ = V3/2. MATLAB simulations are used to predict the 
performance of the prototype antenna. The prototype antenna is shown in Figure 5-1. 





Figure 5 -1: Block Diagram of the Unsealed RSNS Antenna. 
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Every design must begin with one decision. The OSNS system previously 
designed and tested at the Naval Postgraduate School (described in Chapter III) possessed 
two channels and a dynamic range of 66. In order to design a prototype of comparable 
spatial resolution, a two-channel system with dynamic range of 64 was chosen. The 
smaller dynamic range allows the digital hardware to have a 6-bit output, making the 
RSNS system less complex than the OSNS (7-bit) prototype. The second goal of the 
prototype was to demonstrate the ambiguity resolution capability of the Robust 
Symmetrical Number System. 

The dynamic range for a two-channel case is M = 4 m x +lm 2 -2 [14] for moduli 
spaced three or more integers apart (Equation 4.8). Moduli pairs which have a dynamic 
range of 64 are {(1,31), (2,29), (3,27), (4,25), (5,23), (6,21), (7,19), (8,17), (9,15), and 
(10,13)}. The moduli pair (8,17) was chosen to build the prototype. The modulus 8 
channel has two folds and the modulus 17 channel has jy of a fold. This combination of 
folds will nicely demonstrate the RSNS ambiguity resolution capability. Table 5-1 shows 
the vector sequence for mj =8 and m 2 — 17. The dynamic range falls between Vectors 21 
and 84. Note that Vector 20 is ambiguous with Vector 52, and Vector 85 is ambiguous 
with Vector 53. The vectors in the dynamic range are re-indexed with a bin number in 
the 4 th , 8 th , and 12 th columns. 
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Table 5- 1: Moduli 8 and 17 Vector Sequence. 


The selection of moduli determines the required comparator thresholds. As 
determined by Equation 4.12, the comparator thresholds for the two channels are shown 
in Table 5- 2. 
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B. FOLDING WAVEFORMS 

The folding waveforms require two pieces of information: the broadside phase 
difference and the scale factor. Broadside incidence corresponds to the transition 
between Bins 31 and 32. Figure 5-2 shows these bins. In the modulus 8 channel, 
broadside is midway between the comparator #7 and comparator #6 thresholds. Since 
F g 7 = cos(if) and V % 6 = cos(^f), the forced phase difference <f> for the modulus 8 channel 

is = if = f • The forced phase difference for the modulus 17 channel is simpler since 
broadside occurs at the threshold for the ninth comparator. From Equation 4.12, 
V l79 = cos(^f) = cos(f), therefore, <j> x7 = f . 
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Figure 5-2: Vectors around Broadside. 


The equations of the folding waveforms with the correct broadside phase 
difference are (from 4.13): 

mx g (9) = cos(M 8 sin 9 + f), (5.1) 

mx X7 (9) = cos(M 17 sin 9 + f). (5.2) 

The folding waveforms for the prototype antenna (Equations 5.1 and 5.2) are 
simulated for the unsealed array by MATLAB codes rsns8.m and rsnsU.m presented in 
Appendix A. The simulated phase detector outputs are shown as Figure 5-3. The 
BinPop.m algorithm simulates the comparator thermometer code. Its output is shown in 
Figure 5-4. The folding waveforms and thermometer codes for the scaled array are 
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Normalized Mixer Output 


simulated by the same routines and are presented in Figure 5 - 5 and Figure 5 - 6. The 
folding waveforms are frequency independent. 



Angle of Arrival (degrees) 


Figure 5-3: Folding Waveforms for Unsealed Prototype Array. 
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Figure 5-6: Comparator States for the Scaled Prototype Array. 


C. ELEMENT SPACING 

The number of channels, dynamic range, modulus, scale factor, and wavelength 
determine element spacing (Equation 4.20). The first four variables are set. Only the 
wavelength remains variable. The OSNS prototype employed coaxial waveguide 
transitions to receive the incident signal. These components are 1.625 inches in width; 
therefore, the smallest element separation must be greater than 1.625 inches while 
maintaining a frequency in the bandwidth of the component (8.2 to 12.4 GHz). The 
largest modulus determines the minimum element spacing because the modulus is in the 
denominator of Equation (4.21). The calculated element spacings are too small even at 
the longest wavelength: dn(8.2 GHz) = 0.68 inches. To accommodate these components, 
the scale factor would have to be 2.39 or greater. This scale factor corresponds to a 
maximum mappable aperture of 25°. This mappable aperture is unacceptably small. 
Stripline antenna elements were designed and manufactured. These stripline antenna 
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elements are 50 mils thick. The resonant frequency of the stripline antenna elements is 
8.0 GHz. The wavelength is 1.48 inches. The original array did not scale the element 
spacings; however, mutual coupling concerns prompted the construction of a scaled 

array. Using equation 4.20 and £ = V3 /2, the spacing for each element is presented in 
Table 5-3. 


Unsealed Antenna 
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Table 5- 3: Element Spacing for the Prototype Antenna Arrays. 


D. TRANSFER FUNCTION 

The MATLAB algorithm BinPop.m that simulated the thermometer code also 
simulates the transfer function of the entire system. The simulation divides the entire 
180° field of view into 0.1° units. This resolution is the minimum step size of the 
antenna positioning servomotor, so simulated and measured data will be convenient to 
compare. In each unit, the normalized mixer output is calculated for both channels, and 
the number of active comparators is determined by comparing the mixer output to the 
comparator thresholds. The number of active comparators either maps into a bin within 
the dynamic range (0-63) or is assigned an ‘unmappable’ value of 100. The bin mapping 
is shown in Figure 5-7 and Figure 5-8. On the scaled antenna, angles of arrival greater 
than 60°, the RSNS system does not map accurately, as expected. 
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Angle of Arrival (degrees) 

Figure 5 - 7: Unsealed Antenna Map of Angle of Arrival to Bin Number. 
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Figure 5-8: Scaled Antenna Map of Angle of Arrival to Bin Number. 
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The bin number corresponds to an angle of arrival, so the same MATLAB 

A 

algorithm assigns the angle at the bin center 9- to any unit mapping into that bin. For 

angles that map into bin 100, the error bin, an angle of arrival of zero is coded as zero. 
The quantized angle of arrival is the final output of the antenna. Figure 5-9 and Figure 5 
-10 show the transfer functions from angle of arrival to reported angle of arrival. 
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Figure 5-9: Unsealed Antenna Transfer Function. 
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Figure 5 -10: Scaled Antenna Transfer Function. 

1. Bin Population 

Each 0.1° is a population point. The “Bin Population” is the quantity of sample 
points in each bin. BinPop.m counts the bin population as the transfer function is applied. 
Figure 5 - 11 and Figure 5-12 shows the bin population for all 64 bins. A small bin 
population corresponds to a small bin width. This corresponds to a sharp resolution. If 
the system possessed an infinite number of bins, the bin width would be infinitesimally 
small, and a smooth curve would result. The bin population curve resembles a cosecant 
function. The jaggedness of the bin population curves is due to the quantization of the 
population points. The curve becomes smoother as more points are measured and 
fractional quantization error becomes smaller. Finally, the bin population curve of the 
scaled array is only about 85% of the unsealed values. This occurs because the scale 
factor (approximately 0.866) compresses all of the bin widths by a ratio of the scale 
factor. This means the bin is a narrower target for the thermometer code. 
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Figure 5-11: Bin Population of the Unsealed Antenna 
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Figure 5 -12: Bin Population of the Scaled Antenna. 
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Predicted Errors 


The true angle of arrival is a continuous function, but all antenna must map the 
angle of arrival into a finite resolution. This quantization requires that the angle of arrival 
be estimated. High resolution corresponds to small estimation errors. Figure 5-13 and 
Figure 5-14 show the errors that arise due to quantization error. These figures illustrate 
that the antenna reports the angle of arrival to within one degree throughout most of the 
mappable field of view. It also shows that the scaled antenna has smaller estimation 
errors. This confirms the fact that the scaled antenna has a finer resolution and is 
corollary to the smaller bin population. 
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Figure 5-13: Quantization Error in Estimated Angle of Arrival for Unsealed Antenna. 
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Figure 5 -14: Quantization Error in Estimated Angle of Arrival for Scaled Antenna 


E. MUTUAL COUPLING 

The RSNS architecture, like all phase sampling systems, measures the output of a 
phase detector. Phase detectors are efficient in the absence of mutual coupling but 
become inefficient when mutual coupling is present. Practical antenna systems have 
significant levels of mutual coupling. The stripline antenna elements used in the RSNS 
prototype are closely spaced dipoles, so significant mutual coupling is anticipated. 
Mutual coupling causes the amplitude of the signal to vary between elements. This 
amplitude variation is a function of the direction of arrival of the signal. Thus, in the 
presence of mutual coupling, signal amplitude contains direction information. Discarding 
the amplitude information will decrease estimation accuracy [4]. 

The RSNS system operates low noise amplifiers (LNA) in saturation to equalize 
the amplitude imbalances received by each element. The second effect of mutual 
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coupling is to introduce phase errors into the signal due to the presence of secondary 
fields generated by the receiving elements. Although these secondary fields are weak, 
the elements are in close proximity, so the phase error can be significant. The secondary 
fields are created because the receiving elements reflect a portion of the incident wave 
and re-radiate a fraction of the signal due to reflections internal to the microwave circuit. 
The secondary waves are received at the other elements due to the inefficiency of each 
element. Dipoles have simple antenna patterns and the strength of the secondary fields 
can be simulated. The secondary fields also induce reactive fields in the other elements 
that cause another reaction in the other elements. The field interactions are interlaced. 
This is mutual coupling. An antenna simulation can predict the phase error generated by 
mutual coupling. The angle of arrival dependent result is shown as Figure 5-15. 
Element 1 is the reference element, element 2 is the modulus 17 element, and element 3 
is the modulus 8 element. 

Arrays with many elements have lower mutual coupling for the elements near the 
center of the array because the induced fields of the outer elements tend to cancel near the 
array center. Figure 5-16 illustrates an array with an additional pair of elements added 
to the ends of the array. The extra elements do not have circuits attached. Instead, the 
the elements are loaded with a 50Q dummy load that substitutes for the typical 
impedance of microwave circuits. The dummy element spacings are equal to the spacing 
of the adjacent element. The mutual coupling simulation is run for the array with dummy 
elements, and the results are shown in Figure 5-17. The dummy elements reduce the 
mutual coupling induced phase error by approximately 50%. The scaled antenna will be 
designed to incorporate such dummy elements. 
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Figure 5-15: Predicted Mutual Coupling Phase Error for the Unsealed Antenna. 



Figure 5 -16: Unsealed RSNS Array with One Pair of Dummy Elements. 




















Figure 5-17: Simulated Mutual Coupling for Unsealed RSNS Array with Dummy Elements. 

These results predicted in this chapter are the baseline for the prototype arrays. In 
the next chapter the hardware design that implements the RSNS technique will be 
presented. 
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VL RSNS HARDWARE DESIGN 


A. ANALOG CIRCUIT COMPONENTS 

The microwave circuit shown in Figure 6-1 receives the incident waveform using 
stripline antenna elements. The entire microwave circuit is mounted on the back of the 
antenna ground plane. These elements efficiently couple the microwave signal from an 
atmospheric medium into rigid coaxial cables. The signal is amplified using Low Noise 
Amplifiers (LNA) to boost the signal and to provide a signal of constant amplitude for all 
AOA. The low pass filters (LPF) remove all harmonic output from the amplifiers, and 
the bandpass filters remove the DC component of the mixer output. A power splitter 
separates the reference signal into two equal signals to be mixed with the modulus 
signals. The modulus signals are phase shifted to set the phase of the signals at broadside 
relative to the reference. Attenuators are used to match the power from each signal prior 
to mixing. The mixer multiplies the reference and modulus signals together, and an 
internal low pass filter outputs only the low frequency term that is proportional to the 
phase difference between the modulus and reference signals. This signal has a dynamic 
range of approximately ±400 mV. It is provided to the digital circuit (located outside of 
the chamber) via flexible coaxial cables. 

1. STRIPLINE ANTENNA ELEMENTS 

For the moduli and frequency chosen, standard coaxial transition receiving 
elements, such as those used in the OSNS prototype, are too wide for the required 
element spacing. Stripline antenna elements are designed and used to receive the incident 
waves. These elements are shown in Figure 6-2. Their 60 mil width easily fits within the 
element spacing. The printed dipoles are designed to resonate at 8.0 GHz. The 
conceptual linear array is shown in Figure 6- 3. Single dipoles have a uniform receiving 
pattern; however, this pattern is perturbed when the element is located within an array. 
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The array used to measure the element patterns is shown as Figure 6- 4, and Figure 6- 5 
shows the scaled array in the anechoic chamber ready for test. Figure 6- 6 shows the 
measured and predicted element patterns in the array. The measured patterns closely 
match the predicted patterns. The power incident upon each element at broadside is 
approximately -40dBm. 



Figure 6-1: Analog Circuit Diagram Showing Scaled Element Spacings. 
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Figure 6- 2: Stripline Element Construction. 



Figure 6- 3: Stripline Element Array Isometric. 





















































































Figure 6- 5: Scaled Array Mounted on Antenna Pedestal. 
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Figure 6- 6: Stripline Antenna Element Patterns. 


2. Low Noise Amplifiers 

The element patterns in Figure 6- 6 do not have a uniform output over the entire 
AOA. The phase detector requires equal amplitude inputs to operate properly. The low 
noise amplifiers (LNA) provide a uniform amplitude input to the phase detector for all 
AOA. Two LNA are in series after each stripline element. The first amplifier amplifies 
the signal allowing the second amplifier to operate in saturation. The second amplifier 
remains in saturation for received powers as low as -59 dBm as shown in Table 6-1. The 
cascaded amplifier pairs are labeled as shown in Figure 6- 7 
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Line 

Amp 

Make 

Model 

Serial 

Saturated 
Power Out [dBm] 

Minimum Input 
Power to Saturate 

Ref 

R-1 

R-2 

Avantek 

Avantek 

SWL87-1000 

SWL87-1000 

F00553-8728 
F01342-8737 

22.82 

-59 

Mod 8 

8-1 

8-2 

DBS 

DBS 

DWT 18636 
SA84-1043 

006-9923 
F03105-8823 

20.83 

-57 

Mod 17 

17-1 

17-2 

DBS 

Avantek 

DWT 18636 
DWT 18636 

001-9733 

005-9922 

21.74 

t 

Ol 

Oi 


Table 6 - 1 : Cascaded Amplifier Pairs with Saturation Data. 


Reference Modulus 17 Modulus 8 



Figure 6- 7: Cascaded Amplifier Pairs. 


3. Filters 

The amplifiers, operating non-linearly in saturation, produce harmonics as well as 
the carrier frequency. These harmonic signals have power 20dB lower than the 
fundamental frequency. The LNAs also have a DC output. These sidebands and DC bias 
are significant enough to cause mixer phase output errors as large as 20 degrees. Two 
filters are used to suppress these signals. A low pass filter with cutoff frequency of 9 
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GHz is used to attenuate the harmonics. The low pass filter is a MICROLAB/FXR 
Model LA-90N. The stop band attenuation is 55dB minimum. The low pass filter 
outputs were measured. It attenuated the amplifier harmonic signals to less than 90dB 
below the fundamental. A band pass filter with pass band between 6 and 9 GHz is used 
to eliminate the DC bias of the amplifiers. The band pass filter is an AIRTRON part 
number 456969 and has a pass band of 4-8 GHz. These filters effectively eliminate all 
signals except for the fundamental and prevent phase errors at the mixer. 

4. Phase Shifters 

The phase difference between the reference signal and the modulus signal at 
broadside defines the relationship of the folding waveforms. Phase shifters are used to 
“phase in” the antenna. The phase shifter is a variable length coaxial cable. At 8.0 GHz, 
one turn is approximately 8° of phase. Shortening makes the phase more positive, 
lengthening makes the phase more negative. The total range of the phase shifter is 
approximately 80°. 

5. Attenuators 

As shown in Table 6-1, the amplifier pairs have different power outputs. The 
microwave circuitry with all of the components and connectors attenuate the signals 
differently in each channel. To provide signals of equal amplitude into the phase 
detector, the signal amplitudes are measured at the input to the mixer for each channel. 
Attenuators placed at the mixer input equalize the channel inputs. For the prototype, a 
6dB attenuator is placed on the reference line before the modulus 8 mixer and a 2.5dB 
attenuator is placed on the reference line before the modulus 17 mixer. 

6. Phase Detector 

The phase detector is the final component in the microwave circuit. The 
prototype antenna uses two MITEQ model ZMX-10G phase detectors. The phase 
detector is comprised of two parts. The mixer is used to achieve frequency conversion. 
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The information contained in the input signal does not change, only the frequency is 
shifted. Many intermodulation product and harmonic frequencies other than the desired 
frequency exist at the output of a mixer. The high frequency terms (the term giving the 
sum of the input frequencies) are filtered through a low pass filter. The only term, which 
survives out of the low pass filter, is the DC term. 


B. DIGITAL CIRCUIT COMPONENTS 

The digital circuit (Figure 6- 8) maps the two phase detector outputs to an angle 
of arrival. The input from each phase detector varies from approximately -400 mV for 
180° out of phase to +400 mV for signals that are in phase. The signals are shifted and 
amplified in the range 1 to 13 volts. Because the phase signal is always positive, only 
one power supply is required for the reference voltage. The 12 volt difference between 
the high and low values is a large band for the comparator thresholds and uses the power 
supply already present for the comparators. The shifted and amplified phase signal is 
input to all 25 comparators. This signal is compared to the comparator thresholds. Each 
comparator’s TTL output is input to a pin on the Erasable Electronically Programmable 
Read Only Memory (EEPROM) chips. The two EEPROMs map the active comparators 
from each channel to a bin. The six-bit output is sent to the computer via the DIO board. 
The circuit is pictured on top of the computer in the anechoic laboratory in Figure 6- 9 
and Figure 6-10. 
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Figure 6- 9: RSNS Modulus 17 Digital Circuit. 
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Figure 6-10: RSNS Modulus 8 Digital Circuit. 
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Shift and Bias Circuit 


The phase detector output is bipolar and unsuitable for threshold signal 
processing. Its small magnitude (approximately 400mV) also requires a tight reference 
voltage tolerance. The AC component requires the comparators to have multiple power 
supplies for the reference voltage generation. Both of these are undesirable. The shift 
and bias circuit converts the 400 millivolt AC signal into a positive voltage in the range 
of 1-13 volts. This allows the reference voltages to use a significant portion of the 
existing 18V range. The comparator thresholds are spread apart, and only one power 
supply is required to generate the reference voltage. The amplification occurs in two 
steps. First, the phase detector output is amplified using a weighted summing amplifier. 
Second, the signal is inverted. 

The summing amplifier is shown in Figure 6-11. The amplifier is a National 
Instruments LM-741CN Operational Amplifier. This amplifier has two inputs: the mixer 
output (Mod Signal) and V 2 that acts as a bias voltage. R 2 determines the amount of bias 
V 2 provides. R; determines the influence generated by the mixer input. R 3 sets the gain. 
The equation that governs the amplifier operation is the summing amplifier gain equation 
[15]: 



(V v} 

*i | v _i 

< R ] R 2 j 


( 6 . 1 ) 


Because the +18 volt source is available, V 2 is 18 volts. Because one other variable must 
be chosen, R 3 is set to lOkQ. The two remaining resistance values are determined by 
matching the phase detector output maximum and minimum values to the desired 
outputs. The negative sign on the right side of Equation (6.1) indicates that this is an 
inverting amplifier. The phase detector output has no DC bias. A phase detector output 
of zero maps to the middle of the shift and bias band: -7 Volts. When V/= 0, (6.1) reduces 
to 
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( 6 . 2 ) 


R-> — —Rt 


= -10kQ 


18V 

-IV 


= 25.7kQ. 


The value of Rj is determined by the maximum value of the mixer output. For the 
modulus 17 channel, the maximum mixer output is measured as 365 millivolts and the 
output is —13 volts. Ri is determined as follows: 


-K 

-0365V 

( v ~, + n) 


f- 13V 18V ) 

j_ 

v R 3 R-2 J 

UOAQ ' 25.7 kQj 


-0365V 
-600 iiA 


= 609Q. 


(6.3) 



Figure 6-11: Shift and Bias Weighted Summing Amplifier. 
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The maximum mixer output of the modulus 8 mixer is 0.397 millivolts. A similar 
calculation determines Rj=662Q. R 2 is the same in both channels. R] is a lkQ variable 
resistor and R 2 is a variable lOkQ resistor in series with a 22k Q. resistor. 



Pin-Out 

Figure 6- 12: Shift and Bias Inverting Amplifier. 


The output of the summing amplifier must be inverted before transmission to the 
comparators. This is performed by the inverting stage shown in Figure 6- 12. Because 
the magnitude of the summing amplifier is correct, no amplification is required in the 
inverting stage. The equation that governs the output of an inverting amplifier is 



(6.4) 
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To obtain unity gain with invertion, Ri=R 2 . 20kQ resistors are used because precision 
resistors were available. The output of the inverter is +1 to +13 volts as desired. This 
signal is provided to the comparators to determine the level of each signal. 

2. Comparators 

The comparator circuits have 8 and 15 comparators in parallel. The circuit for 
each comparator is shown in Figure 6-13. 


+18 V 



Figure 6-13: Comparator Schematic and Pin-out Diagram. 


The shift and bias circuit output is the input to the non-inverting terminal. The 
second input is a reference voltage generated in the voltage divider network shown in 
Figure 6-14. Since the 18-volt power supply is available as the power supply for the 
comparator, it is used as the source for the voltage divider. The thresholds for the 
comparators, as explained in the previous section, range from 1 to 13 volts, so the 
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conversion from the normalized thresholds to the actual thresholds is V act = 6 V norm + 7. 

The normalized and actual threshold values are shown in columns two and three of 
Table 6-2. Applying Ohm’s law between the source and inverting termininal of 
the comparator, Equation 6.4 shows the relationship between/? and R 3 : 


r _ Vthreshold 

2 R ., 


(6.4) 


The comparators (National Semiconductor LM31 IN) have a maximum input bias 
and offset currents of 100 nA. This current is h in Figure 6-14. In order to limit the 
input offset and bias currents to less than 1% of the current in the voltage divider, the 
total resistance of the network must not exceed 1.5M Q. However, the resistance values 
must be large enough such that the current required by the system remains small. For 
these reasons, Rj is chosen to be 10k Q . The total current for the voltage dividers is 
27.5mA as shown in the lower left comer of 

Table 6-2. Kirchhoff s Current Law applied at the node above R 2 yields / 3 = J x -J 2 . 
The only variables that remain are R 2 and R 3 . Ohm’s law across these resistors dictates: 

R 2 +R 3 = (6.5) 

/■» 


As an initial value, R 2 is chosen to be 5% of R 3 . Once initial values are determined, 
standard resistor values are chosen and R 2 is made a variable resistor. The voltage 
divider sub-circuit is complete and provides a precise threshold voltage to the inverting 
terminal of the comparator. 

When the shift and bias signal is larger than the reference voltage, the comparator 
is ON. In this state, no current is accepted at the output, so the 5 volt signal above the 
output resistor has no voltage drop. The 5 volts is impressed on the output. This voltage 
is passed to the input of the EEPROM. If the reference voltage exceeds the shift and bias 
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signal, the comparator is OFF. In this state, the output sinks current, and the 5 volts are 
dropped across the output resistor. The output of the comparator is shorted to ground 
through the opamp and passed to the input of the EEPROM. 



Figure 6-14: Voltage Divider Circuit. 
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5K 

833 

-0.19509 
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1.000E+04 

842 

2.281 E+02 

2K 

843 
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851 
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772 

1.748E+02 

10K 

773 

3.496E+03 

33K 

-0.18375 

5.89750 


781 

1.000E+04 

782 

2.321 E+02 

5K 

783 

4.641 E+03 

47K 

0.00000 

7.00000 


791 

1.000E+04 

792 

3.031 E+02 

10K 

793 

6.061 E+03 

56K 

0.18375 

8.10250 


7A1 

1.000E+04 

7A2 

3.899E+02 

10K 

7A3 

7.797E+03 

150K//150K 

0.36124 

9.16745 


7B1 

1.000E+04 

7B2 

4.943E+02 

5K 

7B3 

9.886E+03 

100K 

0.52643 

10.15859 


7C1 

1.000E+04 

7C2 

6.170E+02 

20K 

7C3 

1.234E+04 

120K 

0.67370 

11.04217 


7D1 

1.000E+04 

7D2 

7.558E+02 

20K 

7D3 

1.512E+04 

150K 

0.79802 

11.78810 


7E1 

1.000E+04 

7E2 

9.038E+02 

20K 

7E3 

1.808E+04 

180K 

0.89516 

12.37098 


7F1 

1.000E+04 

7F2 

1.047E+03 

50K 

7F3 

2.093E+04 

180K 

0.96183 

12.77095 


7G1 

1.000E+04 

7G2 

1.163E+03 

50K 

7G3 

2.326E+04 

220K 

0.99573 

12.97441 


7H1 

1.000E+04 

7H2 

1.230E+03 

50K 

7H3 

2.459E+04 

220K 



Current Values Input Bias Current 1.00E-07 



ii 

12 

13__ 

811 

1.688E-03 

812 

1.00E-07 

813 

1.688E-03 

821 

1.599E-03 

822 

1.00E-07 

823 

1.599E-03 

831 

1.433E-03 

832 

1.00E-07 

833 

1.433E-03 

841 

1.217E-03 

842 

1.00E-07 

843 

1.217E-03 

851 

9.829E-04 

852 

1.00E-07 

853 

9.828E-04 

861 

7.667E-04 

862 

1.00E-07 

863 

7.666E-04 

871 

6.01 IE-04 

872 

1.00E-07 

873 

6.010E-04 

881 

5.115E-04 

882 

1.00E-07 

883 

5.114E-04 




711 

1.697E-03 

712 

1.00E-07 

713 

1.697E-03 

721 

1.677E-03 

722 

1.00E-07 

723 

1.677E-03 

731 

1.637E-03 

732 

1.00E-07 

733 

1.637E-03 

741 

1.579E-03 

742 

1.00E-07 

743 

1.579E-03 

751 

1.504E-03 

752 

1.00E-07 

753 

1.504E-03 

761 

1.416E-03 

762 

1.00E-07 

763 

1.416E-03 

771 

1.317E-03 

772 

1.00E-07 

773 

1.317E-03 

781 

1.210E-03 

782 

1.00E-07 

783 

1.210E-03 

791 

1.100E-03 

792 

1.00E-07 

793 

1.100E-03 

7A1 

9.898E-04 

7A2 

1.00E-07 

7A3 

9.897E-04 

7B1 

8.833E-04 

7B2 

1.00E-07 

7B3 

8.832E-04 

7C1 

7.841 E-04 

7C2 

1.00E-07 

7C3 

7.840E-04 

7D1 

6.958E-04 

7D2 

1.00E-07 

7D3 

6.957E-04 

7E1 

6.212E-04 

7E2 

1.00E-07 

7E3 

6.21 IE-04 

7F1 

5.629E-04 

7F2 

1.00E-07 

7F3 

5.628E-04 

7G1 

5.229E-04 

7G2 

1.00E-07 

7G3 

5.228E-04 

7H1 

5.026E-04 

7H2 

1.00E-07 

7H3 

5.025E-04 



Table 6-2: Voltage Divider Calculations. 
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3. 


EEPROM 


The 25 comparator states must be processed into a 6-bit angle of arrival. The 
electronically erasable programmable read only memory (EEPROM) performs this 
mapping function. ATMEL Model AT28C256 is used. The EEPROM is a small logic 
block. Each EEPROM converts 15 input bits to an 8-bit output via its internal 
programming. Because of the limit of 15 inputs, two EEPROM must be used to process 
all 25 comparators. The “Thermometer Code” programmed into the EEPROM is shown 
in Table 6-3. The thermometer code is the string instructions to the EEPROM that maps 
input to output. With two EEPROM the thermometer code has two stages as shown on 
the circuit diagrams of Figure 6-15 and Figure 6-16. 

The first stage encodes the 15 least significant (lowest threshold) comparators of 
the modulus 17 channel into a 4-bit representation of the number of active comparators. 
The encoding is a simple binary conversion. The second column of Table 6 - 3 is the 
number of comparators active in the modulus 17 channel in each bin. The fourth column 
is a binary representation of the state of the input pins for the first EEPROM and the fifth 
column is the hexadecimal representation of that binary number. The sixth column is the 
decimal output desired from EEPROM #1 for the given input; it is the decimal number of 
number of active comparators. The seventh column is the hexadecimal representation of 
that digital number. The result is that the number of active comparators is converted to a 
4-bit input. The shift and bias circuit and comparators for Modulus 17 and EEPROM #1 
are physically located on the same circuit board. The output of EEPROM #1 and 
comparators 16 and 17 must be passed to EEPROM #2 on the second circuit board. 
Figure 6-15 illustrates how this signal is passed to EEPROM #2. 
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RSNS Encoding Scheme 


Channel 

17 


10 

10 

1 L 

:vffv 

W&M 

13 

13 

14 _ 

16 

16 

1 7 

-tnM 

HMM 

15 

15 

14 _ 

>Z'i& m. 

Jj&ir 

12 

12 

11 

tiiiiS 

9 

9 

8 

’ -T 

^ 7?P 

6 

6 

5 

STs V 

- =4 ™ 

i'Mini 

3 

3 

2 

£SS2 

. 0 . 

0 

1 

3 

3 

4 

• 

ij 

6 
6 
7 


XICOR #1:10 LSBofMod 17 Channel 


Hex Input Output 


XICOR #2:7 MSB of Modi 7 Channel 


Hex Out 
XICOR #1 



000 0000 0000 0111 
000 0000 0000 0111 
000 0000 0000 1111 


000 000000004111 

000 00000001 lift 

000 OpOO 0^014111 

boo oboo 06111111 
000 0000 0011 1111 
000 0000 0111 1111 

wmmmmmwm 


Input 

XICOR#2 


000 1111 1100 1010 
000 0111 1100 1010 
000 0111 1100 1011 

QQoooWiim ibii 
000 001111001100 


Hex Input 
XICOR#2 



000 0000 0001 1111 
000 0000 0101 1111 
000 0000 0111 1111 

*166060011011111 

;;idS6bSMiji4ii>;: 


000 0001 1100 1111 
000 0011 1100 1111 
000 0011 1100 1110 

iooo^iip^pSWC 

0660111 itbonoi 

0001111 1100.1401 


000 1111 1100 1100 
001 1111 1100 1100 
001 1111 1100 1011 

3 E& mim: 

11061010 

obi 111111001001 
000 1111 1100 1001 
000 1111 1100 1000 

^QQOoin’imioob 
bobbiiinoboii® 
floob.ooii iioooiife! 



Output 

XICOR#2 


0 

1 

2 

PISS-; 

iWMiB 

6 

7 

8 

12 

13 

14 

r " 45 "^ 

18 

19 

20 

24 

25 

26 

: Bzt? ? 

gSfe •> 

30 

31 

32 

T<\34^f 

: r. 


Hex Out 
XICOR#2 





0000011 1100 0110 

03 C 6 

36 

24 

000 0001 1100 0110 

01C6 

37 

25 

000 0001 1100 0101 

01 C 5 

38 

26 


00 C 5 


% ty&t 

Vibob 066041606106 

0004 



, {WO POOOOIOO OIOO : j 

0044 

;r~ 41 Cr 


000 0000 0100 0011 

0043 

42 

2 A 

000 0000 0000 0011 

0003 

43 

2 B 

000 0000 0000 0010 

0002 

44 

2 C 

' ;6 ooq6o6.oioo6qi”6 : ;| 

0042 

^ 45 H 

2 D 

^booobbbbiwoobi 



$ 2 E 

obbobooiiojoboi 1 

usmmsm 


: 2F 

000 oboo 1100 0000 

00 C 0 

48 

30 

000 00011100 0000 

01 CO 

49 

31 

000000111000001 

01 Cl 

50 

32 


03 C 1 

■ : -r'-Si-k ' \ 

~~ y ] 33 - 

000 001 r 11Q0 0010 

- 03 C 2 ; 


34 



I- '1&'£■? 

, 35 

bobbin 1100 obii 

07 C 3 

54 

36 

000111111000011 

0 FC 3 

55 

37 

000111111000100 

0 FC 4 

56 

38 

4 :i@i 4 dQlio 6 ;“" 

1 FC 4 


W 39 

ftii 11000101 

1 FC 5 

'/' ’ ; 58 ; ^ 

m-mu 



£ : . ‘ 5 pcru; 

3 B 

01111111100 0110 

3 FC 6 

60 

3 C 

001111111000110 

1 FC 6 

61 

3 D 

001111111000111 

1 FC 7 

wmmm 

62 

3 E 


Table 6-3: Thermometer Code of EEPROM Programming. 
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EEPROM: ATMEL AT28C256 
#1: Mod 17 LSB 


EEPROM: ATMEL AT28C2J6 
#2 Combiow 



Blue Connector Bare Connector 


Figure 6- 15: Schematic Diagram of Mod 17 Comparator Output to EEPROM#2 Input. 

EEPROM #2 takes the 4-bit output of EEPROM #1 into the least significant bit 
input. Modulus 17 comparators 16 and 17 are the next two input bits. The next eight 
input bits are the modulus 8 comparator inputs. Table 6 - 3 column 8 shows the binary 
input to EEPROM #2 for each bin and column 9 is the hexadecimal representation of this 
input. Column 10 is the digital output and column 11 is the hexadecimal output. The 
hexadecimal columns are required because the programming of the EEPROM is 
conducted in hexadecimal representation. The 6-bit output of EEPROM #2 is the input to 
the DIO board. The DIO board is the interface with the computer and has 50 p ins 
Figure 6- 16 diagrams how the signal is passed to the proper pins of the DIO board. 
LAB VIEW samples the input from the DIO board and records the data. . 
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EEPROM: ATMEL AT28C256 
#2: Combiner 


DIO Block 



Figure 6-16: Schematic Diagram of EEPROM #2 to DIO Board Wiring. 


The LAB VIEW program records the binary output to the DIO board. The transfer 
function is completed when a MATLAB routine ( AOAjrnap.m ) converts the bin number 
to an estimated angle of arrival. 
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VII. TEST AND EVALUATION OF PROTOTYPE HARDWARE 


A. ANECHOIC CHAMBER AND TEST PROCEDURE 

The prototype antenna is tested in the Naval Postgraduate School Microwave 
Anechoic Chamber. The major components and communication paths used for testing 
are shown in Figure 7-1. The test system architecture centerpiece is a 486 computer. 
One of two LABVIEW programs controls the test via a General Purpose Interface Bus 
(GPIB), the Serial port, and a PC-DIO-96 board. The local oscillator that produces the 
radiated signal is located on a Hewlett Packard 8510 Vector Network Analyzer (VNA). 
In the anechoic chamber, a standard gain feed horn transmits the microwave signal 
generated by the local oscillator. The antenna is placed on a pedestal (Figure 7-1) 
approximately 12 feet from the transmitting hom with the microwave circuit mounted on 
the back of the array ground plane. A servomotor is controlled through the computer 
Serial port and rotates the pedestal with a best resolution of 0.1 degrees. The microwave 
circuit continuously analyzes the received signal, and at each resolution step the phase 
detector output is sampled. The phase detector output is carried out of the anechoic 
chamber on RG-59 coaxial cables. The digital circuits process the phase detector output. 
The 6-bit output of the digital circuit is transmitted to the computer via the DIO board. 
Two HP 3478A digital multimeters measure either the phase detector output or the shift 
and bias circuit output. The computer samples the multimeters through the GPIB bus. 
The experiment control and recording equipment are shown in Figure 7-2. 

The data collection and processing technique is illustrated in Figure 7-3. Two 
types of tests are conducted depending upon the parameter of interest. The first test 
measures the output of the phase detector, and other test measures the output of the shift 
and bias amplifier and the binary angle of arrival estimation simultaneously. One 
LABVIEW program is written for each test, so the tests cannot be run concurrently. 
These programs are presented in Appendix B. When the phase detector output is 
measured, the procedure follows the left path on Figure 7-3. The phase detector output 
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Figure 7 - 1: Anechoic Chamber Configuration. 


of each channel is connected to a multimeters via an RG-59 coaxial cable. The 
LAB VIEW VI code DF_mm_R2 samples the multimeters at each servomotor step. The 
data is saved in column format. The four columns of data are (1) antenna position, (2) 
channel 1 (modulus 17) multimeter output, (3) antenna position, and (4) channel 2 
(modulus 8) multimeter output. This data is physically transported from the lab on a 
diskette. The data is opened in Microsoft Excel using the “text-tab delimited” format. 
Manual normalization of the data is performed. The minimum and ma ximum values of 
the phase detector folding waveforms are located. A column is inserted between the 
second and third data columns. The normalized phase detector output Tj iinorm is calculated 
the using the equations: 
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V =(^ m ax+^ m in)/2> 

-VvJ/2, 


Vi, 


UiZl 

A t] 


(7.1) 

(7.2) 

(7.3) 


The normalized output for each channel is placed in the new third column (modulus 17) 
and the sixth columns (modulus 8). 



Figure 7 - 2: Anechoic Chamber Control and Recording Equipment. 
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The normalized phase detector output is saved in tab delimited text format with a “.dat” 
extension. This is the standard input format for the MATLAB processing programs 
provided in Appendix A. The measured phase detector output is visually compared to 
simulated data using the rsns8.m and rsnslJ.m programs. Phase_error.m calculates and 
plots the phase error between a simulated and normalized measured phase detector 
output. The normalized phase detector output is passed through a simulation of the 
digital processing in BinJPop.m. If the output of the shift and bias circuit is connected to 
the multimeters rather than the phase detector output, the same collection and processing 
may be conducted on this signal; however, this specific configuration was not used in this 
thesis. 

The output of the digital circuit follows the path on the right of Figure 7-3. 
Either the phase detector or shift and bias circuit output is connected to the multimeters. 
The binary output of the digital circuit is connected to the computer via the DIO board. 
The LAB VIEW VI code DF5_DIO_R2 samples the multimeters and DIO board at each 
servomotor step. The data is again saved in column format and must be transported, and 
opened as before. No normalization is required for this data, but it must be saved in tab 
delimited text format with a “.dat” extension. Analysis of the data is performed by one 
program, AOAjnap.m. This program plots the folding waveforms and the bin mapping. 
It converts the bin numbers into an estimated angle of arrival and plots the transfer 
function. 

The duration of a test depends upon the desired resolution of the test. For the 
DF_mm_R2 measurement of the analog data, a 0.5° resolution test duration is 25 minutes 
and a high (0.1°) resolution test requires approximately two hours to complete. The low- 
resolution measurements are used for initial phasing of the microwave circuits, and high 
resolution measurements are used for detailed analysis of the antenna performance. The 
digital output sampling program, DF5_DIO_R2, requires 9 hours for a high resolution 
test. The program can conduct low resolution tests, but none were conducted for this 
thesis research. 
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B. UNSCALED ANTENNA RESULTS 


Two antennas have been built and tested. The first antenna is unsealed (£=1). 

The second antenna is scaled for a 60° field of view (£ = V3/2). Both have been tested 
at a resolution of 0.1°. The results from each antenna will be presented individually 
beginning with the unsealed antenna. Figure 7-4 shows the normalized modulus 8 phase 
detector output output. The solid line represents the simulated results and the dashed line 
shows the measured results. The phase adjusters are tuned to achieve the best fit between 
the measured and simulated results. Figure 7-5 shows the same results for the 
normalized modulus 17 phase detector output. 

Figure 7-6 plots the phase difference between the simulated and measured 
results. Phase error, m finds the difference between the arccosines of the simulated and 
normalized phase detector output data of Figure 7 - 4 and Figure 7-5. The comparator 
thresholds are spaced uniformly in phase, so phase errors larger than one comparator step 
in the measured mixer output disrupts the thermometer code. Minimizing the phase error 
is critical to the success of the RSNS architecture. The normalized measured phase 
detector output is processed by the simulated transfer function using BinPop.m in Figure 
7-7. This figure clearly shows that the large phase errors observed in Figure 7-6 
prevent the RSNS system from properly estimating the angle of arrival. The RSNS 
circuit maps the angles of arrival with small phase (-30° < 0 < 0°) errors without 
difficulty. 

The phase detector output is the input of the shift and bias circuit. The variable 
resistors Ri and R 2 in the shift and bias circuit are used to tune the shift and bias circuit. 
The output of the shift and bias circuit is shown as Figure 7-8. The comparator 
networks process the shift and bias output and are subsequently processed by the 
EEPROM. The digital output of the EEPROM #2 is shown as Figure 7-9. This output 
is input to the computer and is processed by AOAjnap.m to determine the estimated 
angle of arrival with the results shown in Figure 7-10. The unsealed antenna maps from 
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approximately 6 = -70° to 6= +65° with significant gaps in the processing between 0° < 
0 < 10° and 35° < 0 < 50°. 



Figure 7-4: Unsealed Antenna, Modulus 8 Normalized Phase Detector Output. 



Figure 7-5: Unsealed Antenna, Modulus 17 Normalized Phase Detector Output. 
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Figure 7-6: Phase Difference Between Simulated and Measured Waveforms, Unsealed Antenna. 



Figure 7-7: Measured Phase Dectector Output through Simulated Digital Circuit. 
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Figure 7-10: Transfer Function of Unsealed Antenna. 
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c. 


SCALED ANTENNA RESULTS 


The phase errors shown in Figure 7 - 6 are attributed in part to mutual coupling as 
predicted in Figure 5-15. The scaled antenna has larger element spacings by a factor of 

l/£ = 2 /a/ 3 and is designed to accommodate dummy elements. The same sequence of 
measurements is conducted with the scaled antenna with no dummy elements installed to 
investigate the effect of scaling only. Figure 7-11 and Figure 7-12 show that the 
folding waveforms are compressed (as expected). The phase error of the array is shown 
in Figure 7-13. The phase error for this array is smaller in magnitude for near broadside 
angles of arrival, but in general, the scaled antenna with no dummy elements has 
approximately the same phase error. The scaling compresses the mappable field of view 

that makes the bin widths smaller by a factor of V3 / 2. This makes the scaled array less 
tolerant of phase errors than the unsealed array. The performance of the measured mixer 
output when passed through the simulated digital circuit is shown in Figure 7-14. The 
shift and bias circuit output shown in Figure 7-15. The digital circuit processes this 
signal with the digital bin mapping result shown in Figure 7-16. The final, processed 
output is the estimated angle of arrival shown in Figure 7-17. As expected, the transfer 
function is steeper. 
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Figure 7 - 12: Scaled Antenna, Modulus 17 Phase Detector Output. 



















of Arrival 



Figure 7-13: Phase Difference Between Simulated and Measured Waveforms, Scaled Antenna. 



Figure 7 - 14: Measured Phase Dectector Output passed through Simulated Digital Circuit 
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Figure 7-15: Shift and Bias Output, Scaled Antenna. 



Figure 7-16: Scaled Antenna Bin Mapping 
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D. DUMMY ELEMENT INFLUENCE ON MUTUAL COUPLING 

The scaled antenna is built to accommodate two pairs of dummy elements that 
repeat the array spacing as shown in Figure 7-18. The dummy elements are loaded with 
50Q loads and are removable. To test the impact of the dummy elements on the mutual 
coupling, the array is tested with no dummy elements installed, with one pair of dummy 
elements installed, and with both pairs of dummy elements installed. The mutual 
coupling-induced phase errors should be reduced by approximately 50% as shown in 
Figure 5-17. The results are shown for each case in Figures 7-20 through 7-22. 



D = Dummy Element 

R = Reference Element 

17 = Modulus 17 Element 

8 = Modulus 8 Element 


Figure 7-18: Scaled Array Construction with Dummy Elements. 
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Figure 7 - 19: Phase Error Results with No Dummy Elements Installed. 



Figure 7 - 20: Phase Error Results with One Pair of Dummy Elements Installed. 
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Figure 7 - 21: Phase Error Results with Two Pair of Dummy Elements Installed. 


Modulus 

No dummy elements 

installed 

One pair 

dummy elements 

Two pair 

of dummy elements 

8 

1.00 

0.86 

0.89 

17 

1.00 

0.93 

0.87 


Table 7-1: Normalized RMS Phase Error for Various Dummy Element Configurations. 


Table 7-1 summarizes the phase errors observed in each modulus for the patterns 
shown in Figures 7-20 through 7-22. The one dummy pair case was expected to show a 
phase error decrease of 50%. The observed decrease was 14% for modulus 8 and 7% for 
modulus 17. This indicates that mutual coupling may not be the dominant source of 
phase error observed in the measured data. 
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E. THREE CHANNEL SCALED ANTENNA RESULTS 


In order to investigate the performance of the RSNS architecture in an 
environment with significantly less mutual coupling, the stripline elements must be 
replaced with channeled waveguide transitions used on the OSNS prototype. These 
components are too wide to fit into the array as designed. Equation 4.21 governs the 
element spacing. The minimum acceptable mappable aperture is 60 degrees, so £ is 

already at its minimum limit (-s/3/2 ). The largest modulus channel defines the minimum 
spacing. In order to increase the spacing, the dynamic range must increase, but from 
(4.11) the dynamic range increases linearly with the modulus in a two-channel case. For 
this reason, a three-channel array was designed and built. The design equations are the 
same as the two-channel array. The geometry, predicted folding waveforms, and 
simulated transfer function are presented in Appendix C. 

The measured phase detector outputs are shown in Figure 7 - 22, Figure 7 - 23, 
and Figure 7 - 24. Phase errors are evident upon a casual visual inspection of the plots. 
These phase errors are too large accurate bin mapping shown in Figure 7-25. The 
observed performance of the antenna in the absence of significant coupling indicates that 
mutual coupling is not the cause of the phase errors observed in the two-channel array. 
The source must be located in the microwave circuit. 
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Figure 7 - 22: Normalized Measured Phase Detector Output, Three Channel Array, Modulus 6. 



Figure 7-23: Normalized Measured Phase Detector Output, Three Channel Array, Modulus 11. 
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Bin Number , Normalized Phase Detector Output 



24: Normalized Measured Phase Detector Output, Three Channel Array, Modulus 17. 



Figure 7-25: Three Channel Array Measured Phase Detector Output Processed by Simulated 

Digital Circuit. 
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The measured results of all three arrays have significant phase errors at the phase 
detector output. Three methods that reduce mutual coupling have been tested: scaling, 
dummy elements, and channeled waveguides. These methods have a small effect on the 
phase errors. Another source of phase error exists within the microwave circuit. Some 
alternative sources of the phase error and methods to reduce these errors are presented in 
the next chapter. 
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VIII. CONCLUSIONS 


A. PROTOTYPE OUTPUT 

The RSNS prototype, like all phase sampling interferometers, measures the phase 
difference between the reference and signal elements. The microwave circuit conditions 
the received signal, and folding waveforms that are a predictable function of the angle of 
arrival are produced by the phase detector. The digital circuit maps the phase difference 
at the mixer output into an estimated angle of arrival. The predicted output of the 
unsealed prototype array is reproduced as Figure 8-1. The measured phase detector 
output contains phase errors between the simulated and predicted results. These phase 
errors cause the thermometer code to incorrectly map the incident wave to the proper 
angle of arrival. The measured transfer function of the unsealed prototype antenna is 
reproduced as Figure 8-2 



Figure 8-1: Simulated Transfer Function of the Unsealed Prototype Antenna. 
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Figure 8-2: Measured Transfer Function of the Unsealed Prototype Antenna. 

Mutual coupling is a common cause of phase errors between closely spaced 
antenna. In order to separate the elements, a scale factor (£) is derived. The scale factor 

is the ratio of the element spacing of the scaled and unsealed arrays. The scale factor 
narrows the mappable field of view of the antenna and increases the resolution within the 

mappable field of view. For a scale factor of £=V3/2, the maximum mappable aperture 
(#a m ) is 60°. The predicted transfer function for the scaled array is reproduced as Figure 

8-3. The mixers output still contains significant phase errors, and the increased 
resolution of the antenna limits the ability to process the errors. The resulting measured 
transfer function is reproduced as Figure 8-4. The steeper slope of the scaled transfer 
function is clearly visible when comparing Figure 8-2 and Figure 8-4. 
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Estimated Angle of Arrival Reported Angle of Arrival (degrees 
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The observed difference in the phase errors between the scaled and unsealed 
antenna does not correspond to the predicted effects of mutual coupling. An array 
comprised of channeled waveguide receiving elements was constructed and tested. These 
elements have significantly less mutual coupling, but the phase errors observed in the 
arrays built with stripline (dipole) elements persisted. This shows that mutual coupling is 
not the primary cause of the observed phase detector phase errors. 

Because the antenna elements are not causing the phase errors, the microwave 
circuit must introduce the errors. Three likely sources are the low noise amplifiers, 
coaxial phase adjusters, and the rigid coaxial cable. The prototype uses three models of 
LNA operating in saturation. - These amplifiers are not phase matched. Variations in the 
amplitude of the incident signal caused by the receiving element patterns may cause 
phase shifts in the amplifier output. Coaxial phase adjusters are a common source of 
noise within microwave circuits. The rigid coaxial cables that transmit the microwave 
signal between components are candidates for three reasons. First, the cable lengths are 
long and subject to thermal expansion and contraction, particularly in the vicinity of the 
Avantek amplifiers that get hotter than the DBS models. Second, mechanical connectors 
are used at the cable junctions. These components often have reflections that introduce 
noise into the circuit. Finally, the cables function are both transmission lines and 
physical support structures for the microwave circuit. The cables are subject to vibration 
on the rotating antenna pedestal and may randomly change the path length within each 
channel. 

The measured results are promising. The RSNS processing is resilient against 
phase errors of about 15 degrees. The mixer output errors are due in part to mutual 
coupling, but the major source is generated within the microwave circuit. Printed 
microwave circuit boards and the use of common components in all channels will reduce 
these sources. When these analog errors are reduced, the RSNS architecture seems to be 
capable of mapping through the remaining mutual coupling induced phase errors. 
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B. FUTURE WORK 


1. Phase Error Reduction 


The sources of phase error in the microwave circuit are controllable. A printed 
circuit board will improve the performance by eliminating the coaxial cabling as an error 
source. When all channels have identical components, phase errors are likely to be 
uniform in each line and will tend to cancel. Using phase-matched amplifiers will 
eliminate the errors caused by using three different amplifier models. The RSNS 
prototype used cascaded comparators to operate in saturation. The implementation of 
constant output amplifiers will reduce the generation of harmonic and intermodulation 
product signals that contribute to the phase error. 

2. Frequency Agility 


The RSNS architecture faces a significant challenge in making the system 
frequency agile. The RSNS system is based upon a fixed broadside phase difference 
between the reference and signal channels. This phase difference is due to the physical 
path difference between the reference and signal channels. For a fixed difference in cable 
lengths, as the frequency varies, the phase difference also varies. Frequency invariant 
phase sampling techniques [1] are being studied. Another approach would use an active 
path length adjusting circuit to maintain the fixed phase difference. 

The scale factor, £» was introduced as the ratio between element spacings for a 
fixed frequency. Because frequency and wavelength are inversely proportional, the scale 
factor is also equal to the reciprocal of the ratio of wavelengths. This insight enables the 
scale factor to be easily used to map angles of arrival for various frequencies incident 
upon an array with fixed element spacings. The scale factor takes the form: 



( 8 . 1 ) 


where/is the design frequency and f is the measured frequency. For example, Figure 8 - 

* 

5, Figure 8-6 and Figure 8-7 show the transfer function for an unsealed antenna for 
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different incident frequencies. Figure 8 - 5 shows the estimated angle of arrival for the 
design frequency of the element spacings. The transfer function uses the entire dynamic 
range. Figure 8 - 6 shows the estimated angle of arrival for a signal with half the design 
frequency. Note that only half of the bins and consequently half of the resolution are 
used. This corresponds to a resolution twice as large as the design frequency. Figure 8 - 
7 shows the transfer function for an incident wave with twice the design frequency. From 
Equation (4.22), the maximum mappable aperture is 0 MM = arcsin(l/£) = arcsin(|) = 30°. 
The estimated angle of arrival is determined from (4.23). 



Figure 8-5: Simulated Transfer Function for Design Frequency. 
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Figure 8-6: Simulated Transfer Function for Half Design Frequency. 



Figure 8-7: Simulated Transfer Function for Twice Design Frequency. 
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3. Multiple Signals 


The prototype can be modified to process signals in a multiple signal environment 
through the use of narrow bandpass filters. The filter must be active to select the target 
signal. Further investigation into the current commercial techniques should allow this 
application to proceed. 


4. Methods To Arrange Elements 

If more than two channels are desired, the chosen modulus may require the 
elements to be spaced close together in order to have the required separation from the 
reference channel. In Figure 8 - 8, a four-channel case is presented. The moduli of 
{4,5,7,9} are chosen. In this case, the modulus 5 spacing would place the element close 
to the modulus 4 element or the modulus 7 element. By mixing the modulus 5 signal 
with the modulus 7 channel instead of the reference channel, the baseline is increased, but 
mutual coupling can be reduced. 



Figure 8-8: Possible Four Channel Geometry. 
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c. 


CONCLUDING REMARKS 


This thesis research has investigated the performance of the Robust Symmetrical 
Number System as a technique to resolve ambiguities in phase sampling interferometry. 
The design equations for comparator thresholds, broadside phase difference, and a scale 
factor are derived. The simulated transfer function demonstrated that the RSNS antenna 
will have a wide instantaneous field of view and high resolution. Mutual coupling that 
adversely effects many phase sampling interferometers was overshadowed by other 
sources of phase error in the RSNS prototypes. The major source of phase errors is 
within the microwave circuit, but these sources can easily be reduced in future prototype 
generations. The simplicity of the RSNS design and absence of moving parts continue to 
make this architecture an attractive means to resolve ambiguities and control errors in a 
phase sampling interferometer. 


Ill 
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APPENDIX A. MATLAB SIMULATION 
AND DATA EVALUATION CODE 


This Appendix presents the codes that simulate the performance of the RSNS 
antenna. The first three codes simulate the antenna performance: RSNS8.m, RSNS17.m, 
and BinPop.m. All three can also be used to compare measured results to the predicted 
results. Two codes specifically designed to analyze the antenna output are also 
presented. These codes are AOAjmap.m and phase_error.m. 


%RSNS8.m WITH SCALING 
% LT Dave Wickersham, 23 March 2000 

% Generates the simulated folding waveform. Two element 
arrays assumed. 

% Can accept normalized mixer output for comparison to 
simulated waveform. 

•s m = modulus; d = spacing (m) ; f = frequency 
clear 
N=2; 

f=8.0e9; C=3e8; wavl=c/f; 

k=2*pi/wavl; 

rad=pi/l80; 

sf=sqrt(3)/2; % sf is the scale factor 


% ffd is forced phase difference to align signal to 
required zero crossing 
ffd=45; %original is 45! 

% number of moduli in array design and their values 
mm=[8 17]; nm=length(mm); 
dispCmoduli are: '),disp(mm) 

% specify a modulus to plot 
m=8; 

% M is Dynamic Range 
M= 64; 

% nf is number of folds 
nf=M/(2*m*N); 

% determines required element spacing 
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d=nf*wavl/sf/2; 

% start, stop, increment angles 
start=-90; stop=90; inc=0.1; 

N=floor((stop-start)/inc)+1; 

for n=1:N; 

thd(n)=start+(n-1)*inc; thr=thd(n)*rad; 
arg(n)=k*d*sin(thr); 

% mixer output 

mx(n)=cos(arg(n)+(ffd)*(pi/180)); 

end 

load p24.dat 

both8=[thd' mx'] 

save simu8.dat both8 -ascii 


figure(1) 
hold 

plot(thd,mx,'k') ; 

plot(p24(:,4),p24(:,6),'--k') 

hold off 

xlabel('Angle of Arrival'),ylabel('Normalized Mixer 

Output'); 

grid 

%title(['Modulus 8, p24, with ffd = ',num2str(ffd)]) 
legend('Simulated','Measured',0) 
axis([-90,90,-1,1]) 
orient landscape 
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%RSNS17 WITH SCALING 


% LT Dave Wickersham, 23 March 2000 

% Generates the simulated folding waveform. Two element 
arrays assumed. 

% Can accept normalized mixer output for comparison to 
simulated waveform. 

% m = modulus; d = spacing (m); f = frequency 

clear 
N=2; 

f=8e9; c=3e8; wavl=c/f; 
sf=sqrt(3)/2; 
k=2*pi/wavl; 
rad=pi/180; 

% ffd is forced phase difference to align signal to 
required zero crossing 
ffd=90; %original is 90! 


m=17; 

% M is Dynamic Range 
M=64; 

% nf is number of folds 
nf=M/(2*m*N); 

% calculation of required element spacing 
d=nf*wavl/sf/2; 

% start, stop, increment angles 
start=-90; stop=90; inc=0.1; 

N=floor((stop-start)/inc)+1; 

for n=l:N; 

thd(n)=start+(n-1)*inc; thr=thd(n)*rad; 
arg(n)=k*d*sin(thr); 

% mixer output 

mx(n)=cos(arg(n)+(ffd)*(pi/180)); 


end 

load p4.dat 
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bothl7=[thd' mx'] 

save simul7.dat bothl7 -ascii 

figure(1) 
hold 

plot(thd,mx,'k'); 

% LABVIEW Output has mod 17 mixer output in second column. 

% Normalization technique puts normalized mixer output in 
third column. 

plot(p4(:,1),p4(:,3),'.g') 
hold off 

xlabel('Angle of Arrival (degrees)'),ylabel('Normalized 

Mixer Output'); 

grid 

title(['Modulus 17, p32, with ffd = ',num2str(ffd)]) 
legend('Modulus 17','Modulus 8', 0) 
axis([-90,90,-1,1] ) 
orient landscape 


116 


%BinPop.m 

% LT Dave Wickersham 23March2000 

% step size is 0.1 degrees based on minimum stepper motor 
resolution 

% simulates the entire transfer function 

% 1- determines number of active comparators at each step 

for both modulus 
% 2- determine bin number 

% 3- determine angle of arrival 

% Also calculates bin population and quantization error. 

% Load either the simulated data (simul7.dat, simu8.dat) 

% or the experimental data - SEE BELOW (ALSO SEE LINE 155- 
156) 

% Uses predicted thresholds calculated from cosine of 
equally spaced angles 

^*^** , a*r'* , ** i '* , '* , * > *‘* , * , ***'* , * , ^*r-* , ** , tfc' , *'** , ** , ** , *rfr'* , *** , afc , ^fc'‘*''* , -^r 

format long; 
clear; 

sum_error_squared=0; Excess_RE=0; 
sf=sqrt(3)/2; 

load simul7.dat 
load simu8.dat 
dfl=simul7; 

%load p32.dat 

%df1=[p32(:,1) p32(:,3)]; 

%********************************* 
point_num=1801; 

Sum_AE=0; 

%comparators mod 17 

for i=l:point_num 

c_ml7(i,[1])=i; 

C_ml7 (i,[2])=0; 
c ml7 (i, [3])=0; 
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C_ml7(i, [4])=0; 
c_ml7(i, [5])=0; 
c_ml7(i, [6])=0; 
c_ml7 (i, [7]) =0 ; 
c_ml7(i, [8])=0; 
c_ml7(i, [9])=0; 

C_ml7(i, [10])=0; 

C_ml7(i, [11])=0; 

C_ml7(i, [12])=0; 
c_ml7(i, [13])=0; 

C_ml7(i, [14])=0; 

C_ml7(i, [15])= 0; 
c_ml7(i, [16])=0; 

C_ml7(i, [17])= 0; 
c_ml7(i, [18])=0; 

end 

%SIMULATION THRESHOLDS 
tl = cos (33*pi/34) ; 
t2 = cos(3l*pi/34); 
t3 = cos(29*pi/34); 
t4 = cos(27*pi/34); 
t5 = cos(25*pi/34); 
t6 = cos(23*pi/34); 
t7 = cos(21*pi/34); 
t8 = cos(19*pi/34); 
t9 = cos(17*pi/34); 
tlO = cos(15*pi/34); 
til = cos(13*pi/34); 
tl2 = COS(ll*pi/34); 
tl3 = COS(9*pi/34); 
tl4 = cos (7*pi/34) ; 
tl5 = cos(5*pi/34); 
tl6 = cos(3*pi/34); 
tl7 = COS (pi/34); 

for i=l:point_num 

% all the data points 

% cl is the bottom comparator 

% tl is the threshold voltage for cl etc. 

if (dfl (i, [2]) ) > tl 
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c_ml7(i,[2]) = 1; 
else 

c_ml7 (i, [2] ) = 0; 

end 

if (df 1 (i, [2] ) )> t2 
c_ml7 (i,[3]) = 1; 
else 

C_ml7(i,[3]) = 0; 

end 

if (dfl(i,[2]))> t3 
c_ml7(i,[4]) = 1; 
else 

C_ml7(i,[4]) = 0; 

end 

if (dfl(i,[2]))> t4 
c_ml7(i,[5]) = 1; 
else 

C_ml7(i,[5]) = 0; 

end 

if (dfl (i, [2] ) )> t5 
c_ml7(i,[6]) = 1; 
else 

c_ml7(i,[6]) = 0; 

end 

if (dfl (i, [2]))> t6 
C_ml7(i,[7]) - 1; 
else 

C_ml7(i,[7]) = 0; 

end 

if (dfl(i,[2]))> t7 
c_ml7(i,[8]) = 1; 
else 

c_ml7(i,[8])= 0; 

end 

if (dfl(i,[2]))> t8 
c_ml7(i,[9]) = 1; 
else 

c ml7 (i, [9] ) = 0; 



end 


if (dfl (i, [2] ) ) > 19 
c_ml7(i, [10]) = 1; 
else 

c_ml7 (i, [10] ) = 0; 
end .. 

if (dfl(i, [2]))> tlO 
c_ml7(i, [11] ) = 1; 
else 

c_ml7(i, [11]) = 0; 

end 

if (dfl (i, [2] )) > til 
c_ml7(i, [12]) = 1; 
else 

c_ml7(i, [12]) = 0; 

end 

if (dfl (i, [2] ) ) > tl2 
C_ml7(i, [13]) = 1; 
else 

C_ml7(i, [13]) = 0; 

end 

if (dfl(i,[2]))> tl3 
c_ml7(i, [14]) = 1; 
else 

c_ml7(i, [14]) = 0; 

end 

if (dfl (i, [2]))> 114 
c_ml7(i, [15]) = 1; 
else 

c_ml7(i, [15]) = 0; 

end 

if (dfl(i, [2]))> tl5 
c_ml7 (i, [16] ) = 1; 
else 

c_ml7(i, [16]) = 0; 

end 

if (dfl(i, [2]))> 116 
c_ml7(i, [17]) = 1; 
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else 

c_ml7(i, [17]) = 0; 

end 

if (dfl (i, [2]))> 117 
C_ml7(i, [18]) = 1; 
else 

C_ml7(i, [18]) = 0; 

end 


end 

%encoding the thermometer code into a decimal number. 
%rsns_ml7 contains the number of comparators on at 
each time 

%step 

for j =1:point_num 

A=c_ml7 (j , [2] ) ; 

B=c_ml7 (j , [3] ) ; 

C=c_ml7 (j , [4] ) ; 

D=c_ml7(j, [5]); 

E=c_ml7(j, [6]); 

F=c_ml7(j, [7] ) ; 

G=C_ml7 (j , [8] ) ; 

H=c_ml7(j, [9]); 

I=c_ml7 (j , [10] ) ; 

J=c_ml7 (j , [11] ) ; 

K=c_ml7(j, [12]); 

L=c_ml7(j, [13]) ; 

M=c_ml7(j, [14]) ; 

N=c_ml7(j, [15]); 

0=c_ml7(j, [16]); 

P=C_ml7(j, [17]) ; 

Q=c_ml7(j, [18]) ; 

Z=A+B+C+D+E+F+G+H+I+J+K+L+M+N+0+P+Q; 
rsns_ml7(j)=Z; 

end 

xarray=[-90:180/1800:90]; 
figure(1) 


121 



plot(xarray,rsns_ml7(:)) ,grid 
xlabel('Angle of Arrival (degrees)') 
ylabel('Number of Comparators') 

% title('Number of Comparators on each time step for mod 
17') 

axis([-90 90 -1 17]) 

%********************************* 
clear dfl 

dfl=simu8; 

%df1=[p32(:,4) p32(:, 6) ]; 
%********************************* 

%comparators mod 8 

for i=l:point_num 

c_m8(i,[1])=i; 
c_m8(i,[2])=0; 
c_m8(i, [3])=0; 
c_m8(i, [4])=0; 
c_m8(i, [5])=0; 
c_m8(i,[6])=0; 

C_m8 (i, [7])=0; 
c_m8 (i, [8] )=0; 
c_m8 (i, [9])=0; 

end 

tl = cos(15*pi/l6); 
t2 = cos(13*pi/l6); 
t3 = cos(1l*pi/16); 
t4 = cos(9*pi/l6) ; 
t5 = cos(7*pi/l6) ; 
t6 = cos(5*pi/l6) ; 
t7 = cos(3*pi/l6) ; 
t8 = cos(pi/16); 

for i=l:point_num 

% all the data points 

% cl is the bottom comparator 

% tl is the threshold voltage for cl etc 
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if (df 1 (i, [2] )) > tl 
c_m8(i,[2]) = 1; 
else 

C_m8 (i, [2] ) = 0; 

end 

if (df 1 (i, [2] ) ) > t2 
C_m8 (i, (3] ) = 1; 

else 

c_m8(i,[3]) = 0; 

end 

if (df 1 (i, [2] )) > t3 
C_m8 (i, [4] ) = 1; 

else 

C_m8(i,[4]) = 0; 

end 

if (df 1 (i, [2] ) ) > t4 
c_m8(i,[5]) = 1; 
else 

c_m8 (i, [5] ) = 0; 

end 

if (df 1 (i, [2] ) ) > t5 
c_m8(i,[6]) = 1; 
else 

c_m8 (i, [6] ) = 0 ; 

end 

if (df1(i,[2]))> t6 
C_m8(i,[7]) = 1; 

else 

c_m8 (i, [7] ) = 0; 

end 

if (df1(i,[2]))> t7 
c_m8(i,[8]) = 1; 
else 

c_m8 (i, [8] ) = 0 ; 

end 

if (dfl (i, [2] ) ) > t8 
c m8 (i, [9] ) = 1; 


123 



end 


else 

c_m8(i, [9]) = 0; 

end 


%encoding the thermometer code into a decimal number. 
%rsns_m8 contains the number of comparators on at 
each time 

%step 

for j=l:point_num 

A=c_m8 (j , [2] ) ; 

B=c_m8 (j , [3]) ; 

C=c_m8 (j , [4] ) ; 

D=c_m8 (j , [5]) ; 

E=c_m8 (j , [6] ) ; 

F=c_m8 (j, [7]) ; 

G=c_m8 (j , [8] ) ; 

H=c_m8(j, [9]) ; 

Z=A+B+C+D+E+F+G+H; 

rsns_m8(j)=Z; 


end , 
figure(2) 

plot(xarray,rsns_m8(:),'r'),grid 
xlabel('Angle of Arrival (degrees)') 
ylabel('Number of Comparators') 

% title('Number of Comparators on each time step for m8') 
axis([-90 90 -18]) 

figure(3) 
hold on 

plot(xarray,rsns_ml7(:),'k') 
plot(xarray,rsns_m8(:),'. -k') 
hold off 

xlabel('Angle of Arrival (degrees)') 
ylabel('Number of Comparators') 
legend('Mod 8','Mod 17',0) 

%title('Number of Comparators on each time step for mod 17 
Sc mod 8 ') 

axis([-90 90 0 17] ) 
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orient landscape 


% integers for trans_funct file 

bin=-l; bin_ctr=l; start_th=-90; end_th=0; 

sum_error_squared=0; 

% Bin mapping for thermometer code, doa is the bin number, 
for s=l:point_num 

if (rsns_ml7(s)==10 & rsns_m8(s)==6) 
doa(s) = 0; 

elseif (rsns_ml7(s)==10 & rsns_m8(s)==5) 
doa(s) = 1; 

elseif (rsns_ml7(s)==11 & rsns_m8(s)==5) 
doa(s) = 2; 

elseif (rsns_ml7(s)==11 & rsns_m8(s)==4) 
doa(s) = 3; 

elseif (rsns_ml7(s)==12 & rsns_m8(s)==4) 
doa(s) =4 ; 

elseif (rsns_ml7(s)==12 & rsns_m8(s)==3) 
doa(s) = 5; 

elseif (rsns_ml7(s)==13 & rsns_m8(s)==3) 
doa(s) = 6; 

elseif (rsns_ml7(s)==13 & rsns_m8(s)==2) 
doa(s) = 7; 

elseif (rsns_ml7(s)==14 & rsns_m8(s)==2) 
doa(s) = 8; 

elseif (rsns_ml7(s)==14 & rsns_m8(s)==1) 
doa(s) = 9; 

elseif (rsns_ml7(s)==15 & rsns_m8(s)==1) 
doa(s) = 10; 

elseif (rsns_ml7(s)==15 & rsns_m8(s)==0) 
doa(s) = 11; 
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elseif (rsns_ml7(s)==16 & rsns_m8(s)==0) 
doa(s) = 12; 

elseif (rsns_ml7(s)==16 & rsns_m8(s)==1) 
doa(s) = 13; 

elseif (rsns_ml7(s)==17 & rsns_m8(s)==1) 
doa(s) =14; 

elseif (rsns_ml7(s)==17 & rsns_m8(s)==2) 
doa(s) = 15; 

elseif (rsns_ml7(s)==16 & rsns_m8(s)==2) 
doa(s) = 16; 

elseif (rsns_ml7(s)==16 & rsns_m8(s)==3) 
doa(s) = 17; 

elseif (rsns_ml7(s)==15 & rsns_m8(s)==3) 
doa(s) = 18; 

elseif (rsns_ml7(s)==15 & rsns_m8(s)==4) 
doa(s) = 19; 

elseif (rsns_ml7(s)==14 & rsns_m8(s)==4) 
doa(s) = 20; 

elseif (rsns_ml7(s)==14 & rsns_m8(s)==5) 
doa(s) = 21; 

elseif (rsns_ml7(s)==13 & rsns_m8(s)==5) 
doa(s) =22; 

elseif (rsns_ml7(s)==13 & rsns_m8(s)==6) 
doa(s) = 23; 

elseif (rsns_ml7(s)==12 & rsns_m8(s)==6) 
doa(s) = 24; 

elseif (rsns_ml7(s)==12 & rsns_m8(s)==7) 
doa(s) =25; 

elseif (rsns_ml7(s)==11 & rsns_m8(s)==7) 
doa(s) =26; 


126 



elseif (rsns_ml7(s)==11 
doa(s) = 27; 

elseif (rsns_ml7(s)==10 
doa(s) = 28; 

elseif (rsns_ml7(s)==10 
doa(s) = 2 9; 

elseif (rsns_ml7(s)==9 & 
doa(s) = 30; 

elseif (rsns_ml7(s)==9 & 
doa(s) = 31; 

elseif (rsns_ml7(s)==8 & 
doa(s) = 32; 

elseif (rsns_ml7(s)==8 & 
doa(s) = 33; 

elseif (rsns_ml7(s)==7 & 
doa(s) = 34; 

elseif (rsns_ml7(s)==7 & 
doa(s) = 35; 

elseif (rsns_ml7(s)==6 & 
doa(s) = 36; 

elseif (rsns_ml7(s)==6 & 
doa(s) = 37; 

elseif (rsns_ml7(s)==5 & 
doa(s) = 38; 

elseif (rsns_ml7(s)==5 & 
doa(s) =39; 

elseif (rsns_ml7(s)==4 & 
doa(s) =40; 

elseif (rsns_ml7(s)==4 & 
doa(s) =41; 


& rsns_m8(s)==8) 
Sc rsns_m8 (s) ==8) 
Sc rsns_m8 (s) ==7) 
rsns_m8(s)==7) 
rsns_m8(s)==6) 
rsns_m8(s)==6) 
rsns_m8(s)==5) 
rsns_m8(s)==5) 
rsns_m8(s)==4) 
rsns_m8(s)==4) 
rsns_m8(s)==3) 
rsns_m8(s)==3) 
rsns_m8 (s) ==2) 
rsns_m8 (s) ==2) 
rsns m8(s)==1) 
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elseif (rsns_ml7(s)==3 & rsns_m8(s)==1) 
doa(s) = 42; 

elseif (rsns_ml7(s)==3 & rsns_m8(s)==0) 
doa(s) = 43; 

elseif (rsns_ml7(s)==2 & rsns_m8(s)==0) 
doa(s) = 44; 

elseif (rsns_ml7(s)==2 & rsns_m8(s)==1) 
doa(s) = 45; 

elseif (rsns_ml7(s)==1 & rsns_m8(s)==1) 
doa(s) = 46; 

elseif (rsns_ml7(s)==1 & rsns_m8(s)==2) 
doa(s) = 47; 

elseif (rsns_ml7(s)==0 & rsns_m8(s)==2) 
doa(s) =48; 

elseif (rsns_ml7(s)==0 & rsns_m8(s)==3) 
doa(s) = 49; 

elseif (rsns_ml7(s)==1 & rsns_m8(s)==3) 
doa(s) = 50; 

elseif (rsns_ml7(s)==1 & rsns_m8(s)==4) 
doa(s) = 51; 

elseif (rsns_ml7(s)==2 & rsns_m8(s)==4) 
doa(s) = 52; 

elseif (rsns_ml7(s)==2 & rsns_m8(s)==5) 
doa(s) = 53; 

elseif (rsns_ml7(s)==3 & rsns_m8(s)==5) 
doa(s) = 54; 

elseif (rsns_ml7(s)==3 & rsns_ra8(s)==6) 
doa(s) =55; 

elseif (rsns_ml7(s)==4 & rsns_m8(s)==6) 
doa(s) = 56; 
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elseif (rsns_ml7(s)==4 & rsns_m8(s)==7) 
doa(s) = 57; 

elseif (rsns_ml7(s)==5 & rsns_m8(s)==7) 
doa(s) = 58; 

elseif (rsns_ml7(s)==5 & rsns_m8(s)==8) 
doa(s) =59; 

elseif (rsns_ml7(s)==6 & rsns_m8(s)==8) 
doa(s) = 60; 

elseif (rsns_ml7(s)==6 & rsns_m8(s)==7) 
doa(s) = 61; 

elseif (rsns_ml7(s)==7 & rsns_m8(s)==7) 
doa(s) = 62; 

elseif (rsns_ml7(s)==7 & rsns_m8(s)==6) 
doa(s) = 63; 

else 

doa(s) = 100; % if the vector does not map, 100 

is the error code 

end 


% Uses the formula asin[(#-34)/64] to find bin center 
gamma=2*doa(s)+1; 

bin_center(s)=asin(gamma*sf/64)*180/3.1415; 

% if 

if doa(s)==100 

bin_center(s)=0; % if the vector did not map, 

reports as broadside 
end 

angle_error(s)=bin_center(s) -(s-900)/10; 

%Calculation of RMS Reporting Error 

%if statement eliminates excessive errors and subs with 
mean error. 

if abs(angle_error(s)) <= 15 

Sum_AE = Sum_AE + angle_error(s)*angle_error(s); 
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end 


Exc e s s_RE=Exce s s_RE+1 ; 

end 


z=0:63 ; 

for 1=1:64 

bin_count(1)=0; 

end 

for lb =1:1801 

if doa(lb)==100 
bc=l; 
else 

bc=doa(lb)+1; 

end 

bin_count(be)=bin_count(be)+1; 

end 


figure(4) 
hold 

plot(z,bin_count,'k') , axis([0,63,0,60]),grid 
hold off 

xlabel( 1 Bin Number') 
ylabel('Bin Population') 


figure (5) 

plot(xarray,doa,'k'),axis([-90,90,0,65]) 
xlabel('Angle of Arrival- (degrees)') 
ylabel( 1 Bin Number') 
grid 

figure(6) 

plot(xarray,bin_center,'k'),grid 

axis([-90,90,-90,90]); 

xlabel('Angle of Arrival (degrees)') 

ylabel('Reported Angle of Arrival (degrees)') 


figure(7) 

plot(xarray,angle_error,'k'),grid 

axis([-60,60,-2,2]); 

ylabel('Reporting Error (degrees)') 
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xlabel('Angle of Arrival (degrees)') 
end of program. 
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% AOA_Map.m Maps bin number to the angle of arrival for 
measured data 

% MAKE SURE THE MEASURED BIN IS IN THE 7TH COLUMN! ! ! 

% LT Dave Wickersham 16 Mar 00 


clear all; 

M=64; % dynamic range 

sf=1;%sgrt(3)/2; % scale factor 

load p30.dat; 

dfl=p30; 

doa=df1(:,7) ; 

point_num=length(df1); % determines 

number of samples 

xarray=[-90:180/(point_num-l):90] ; % determines 

resolution bit 


sbl7=df1( :,2) ; 

Sb8=df1(:,4); 

for s=l:point_num 

gamma(s)=-l+(2*doa(s)+1)/M; 
bin_center(s)=asin(gamma(s)/sf)*180/pi; 
if doa(s)==100 

bin_center(s)=0; 

end 

end 

figure(1) 

plot(xarray,doa, 'k'),axis([-90,90,0,65] ) 

% title('Bin Number for p33') 

xlabel('Angle of Arrival (degrees)') 
ylabel('Bin Number') 
grid 

figure(2) 
hold 

plot(xarray,bin_center,'k'),axis([-90,90,-90,90]) 
% title('Estimated Angle of Arrival for p33') 
xlabel('Angle of Arrival (degrees)') 
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ylabel('Estimated Angle of Arrival') 

grid 

hold off 

figure(3) 
hold 

plot(xarray,sb8,'k') 

plot(xarray,sbl7,* —k') 

axis([-90,90,1,13] ) 

legend('Modulus 8','Modulus 17',0) 

hold off 

xlabel('Angle of Arrival (degrees)') 
ylabel('Shift and Bias Output (volts)') 
' grid 


%phase_error.m Determines the phase error between a 
predicted and measured waveform 
%2 channel input 

clear all; 

load p9.dat 
theta= [p9 (:, 1) ] ; 

dv8=[p9 (:,6)] ; 
dvl7= [p9(:,3)] ; 
md=length(dv8); 

% constants and design variables 
N=2 ; 

f=8e9; c=3e8; wavl=c/f; 
sf=2/sqrt(3); 
k=2*pi/wavl; 

M=64; 

nf=M/(2*N);% nf is number of folds (not really, no modulus 
in this equation) 

dw=nf/2; % new design for mixer ouputs 

FirstMin=0;FirstMax=0;SecMin=0;SecMax=0;Maxi7=0;Minl7=0; 

rad=pi/180; 

thet r=theta*rad; 
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ffd8=175; %45 

ffdl7=7 0; %90 

e 8_sq=0; el7_sq=0; cnt8=0; cntl7=0; 
d8=dw*wavl*sf/8; 
dl7=dw*wavl*sf/17 ; 

%subroutine to find local max and min 

% works because I know about where the max and min will 
fall 

% used to get the right quadrant for the arccos 
%Mod 8 

for i=(md-1)/6:(md-1)/2 
if dv8(i)<FirstMin 
FirstMin=dv8(i); 

IndexlMin=i; 

end 

end 

for i=(md-1)/3:(md-1)/2 
if dv8(i)>FirstMax 
FirstMax=dv8(i); 

IndexlMax=i; 

end 

end 

for i=(md-1)/2:2*(md-1)/3 
if dv8(i)<SecMin 
SecMin=dv8(i); 

Index2Min=i; 

end 

end 

for i=2*(md-1)/3:8*(md-1)/9 
if dv8(i)>SecMax 
SecMax=dv8(i); 

Index2Max=i; 

end 

end 

%Mod 17 

for i=(md-1)/2:5*(md-1)/6 
if dvl7(i)<Minl7 
Minl7=dvl7(i); 

Indexl7Min=i; 
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end 


end 

for i=(md-1)/6:(md-1)/2 
if dvl7(i)>Maxl7 
Maxl7=dvl7(i); 
Indexl7Max=i; 

end 

end 


for n=l:md 

delt8(n)=k*d8*sin(thet_r(n))+(ffd8)*rad; %delt8 is the 
simulated total angle in rad 
unwrap(delt8(n)); 

%while delt8(n)>pi 
% delt8(n)=delt8(n)-2*pi; 

%end 

%while delt8(n)<-pi 
% delt8(n)=delt8(n)+2*pi; 

%end 

dm8(n)=acos(dv8(n)); %dm8 is measured data's total 
angle in rad 
% unwrap(dm8(n); 

if (n>=IndexlMin)&(n<=IndexlMax) 
dm8(n)=-dm8(n); 
end 

if (n>=Index2Min)&(n<=Index2Max) 
dm8 (n) =-dm8 (n) ; 
end 


phi8(n)=(delt8(n)-dm8(n)); 
% unwrap phi8(n); 
while phi8(n)>pi 

phi8(n)=phi8(n)-2*pi; 

end 

while phi8(n)<-pi 

phi8(n)=phi8(n)+2*pi; 

end 
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deltl7(n)=k*dl7*sin(thet_r(n))+(ffdl7)*rad; %deltl7 is 
the simulated total angle in rad 

% unwrap(deltl7(n)); 

while deltl7(n)>pi 

deltl7(n)=deltl7(n)- 2 *pi ; 

end 

while deltl7(n)<-pi 

deltl7(n)=deltl7(n)+2*pi; 

end 

dml7(n)=acos(dvl7(n)); %dml7 is measured data's total 
angle in rad 

if n>=Indexl7Min 

dml7(n)=-dml7(n); 

end 

if n<=Indexl7Max 

dml7 (n) =-dml7 (n) ; 

end 


phi17(n)=(deltl7(n)-dml7(n)); 
while phil7(n)>pi 

phi17(n)=phil7(n)-2*pi; 

end 

while phil7(n)<-pi 

phil7(n)=phil7(n)+2*pi; 

end 


end 

phi8 =phi8/rad; 
phi17 =phi17/rad; 


f° r i=l:md %phi contributes to 

RMS error only if it is real 

e8_sg = e8_sq + (real(phi8(i))*real(phi8(i))); 
el7_sq = el7_sq+ (real(phil7(i))*real(phil7(i))); 

end 
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rms8=sqrt(e8_sq)/(md-1) ; %calculate RMS phase error 

for each channel 
rmsl7=sqrt(el7_sq)/(md-1); 

disp([ 1 RMS Phase Error of Mod 8 is ',num2str(rms8),' ']) 
disp(['RMS Phase Error of Mod 17 is ',num2str(rmsl7),' ']) 


figure(1) 
hold 

%plot(pl2(:,4),pl2(:,8), '--r') 
plot(theta,phi8,'k') ; 
plot(theta,phil7,'.-k'); 
axis([-90,90,-90,90] ) 
hold off 

legend('Mod 8','Modl7',0) 

xlabel('Angle of Arrival (degrees) 1 ),ylabel('Phase Error 
(degrees)'); 
grid 

%title(['Phase Error between Simulated and Measured Data, 
Pattern 32 1 ]) 
orient landscape 
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APPENDIX B. LAB VIEW VI CODE 


LAB VIEW Virtual Instruments (VI) code has been used to create the data 
collection programs. These VI codes are updated for the LAB VIEW 5.1 revision and are 
located in the D:\Anechoic Chamber\Programs\vi for LV5.1\RSNS-OSNS Antenna 
directory. 

The General Purpose Interface Bus (GPIB), also called IEEE 488, is a method of 
communicating with standard-alone instruments, such as multimeters and oscilloscopes. 
The most direct method is to install a plug-in GPIB board in the computer and connect 
the respective instrument directly to this board with a GPIB cable, as the Figure B-l 
shows. The GPIB cable consists of a shielded 24-conductor cable, which has both a plug 
and a receptacle connector at each end, as Figure B-2 shows. 



Figure B-l: Configuration to Link Devices Using a GPIB Cable. [16] 
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The GPIB carries device-dependent messages and interface messages. GPIB 
devices can be Talkers, Listeners, and Controllers. A multimeter, like in our case, is a 
Talker and may be a Listener as well. The Controller manages the flow of information on 
the GPIB by sending commands to all devices. 

With the GPIB interface board and its software, the computer plays all three roles: 

• Controller - to manage the GPIB, 

• Talker - to send data, 

• Listener - to receive data. 
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Figure B-2: A GPIB Connector Showing Signal Assignment. [16] 
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After this small analysis, the code will be understandable easily. 

Any VI contains an interactive user interface, which is called the front panel, 
because it simulates the panel of physical instrument. The front panel can contain knobs, 
push buttons, graphs, and other controls and indicators, as shown in Figure B-3. VI 
receives instructions from a block diagram, which it can be constructed using the 
graphical programming language, G. We can think of the block diagram as source code. 

The procedure for this code and the small steps, that are followed for the 
successful communication between a GPIB device and a computer, are shown below: 

i. Service Request (SRQ). The GPIB line that the device asserts to notify the 
CIC (Controller-In-Charge) that the device needs servicing. 

ii. Writes command string and send it to the GPIB device (multimeter), as the 
frame 0 in Figure B.4 shows. 

iii. Waits for GPIB to assert SRQ to indicate that reading is available, as 
shown in Figure B-5 (frame 1). 

iv. Reads measurement from multimeter, as shown in Figure B.5 (frame 2). 
As the left icon in this frame indicates, reads by the count number of bytes 
from the GPIB device (multimeter). As the right icon indicates, it copies 
the measurement value. 

v. Figure B-6 shows that it clears the SRQ mode in order to be ready for the 
next request. [16] 
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Figure B-4: Writes Command String to the GPIB Device (Multimeter). 
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Figure B-5: Frame 1 Waits for GPIB to Assert SRQ. 
Frame 2 Reads Measurement from Multimeter. 
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Figure B-6: Clears the SRQ Mode to be Ready for the Next Request. 
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APPENDIX C. THREE CHANNEL RSNS 
ARRAY DESIGN RESULTS 


This Appendix presents the design of a 3-channel RSNS prototype array. The 
element geometry is presented and followed by the predicted folding waveforms and 
transfer function. The measured results are presented in Chapter VII Section D. 

The moduli are mj = 6 , m 2 = 11, and m 3 = 17. The dynamic range is 226. The 
element spacings are calculated from (4.20) with N- 3. The two dummy elements are 
placed to minimize any mutual coupling that may be present. The array geometry is 
shown in Figure C-l. 



Figure C-l: 3 Channel RSNS Prototype Array Geometry. 


The folding waveforms for each modulus are shown as Figures C-2 through C-4. 
The thermometer code is presented as Table C-l with the comparator thresholds 
determined by (4.12). Figure C-5 shows the simulated bin mapping of the antenna. The 
RMS resolution is 0.55 degrees. 
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Figure C-2: Folding Waveform for Three Channel Array Modulus 6. 



Figure C-3: Folding Waveforms for Three Channel Array Modulus 11. 
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Bin Number Normalized Phase Detector Output 



Figure C-4: Folding Waveforms for Three Channel Array Modulus 17. 



Figure C-5: Bin Mapping of RSNS 3 Channel Prototype Array. 
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